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Research  during  the  reporting  period  was  accomplished 
in  four  categories  —(a)  Analysis  of  Rayleigh  wave 
dispersion  for  polar  paths;  (b)  development  of  improved 
array  processing  techniques;  (c)  determination  of  bias  in 
network  computations  of  body  wave  magnitudes;  and  (d) 
investigation  of  the  effect  on  the  Rayleigh  spectrum  of 
lateral  heterogeneity  in  earthquake  source  regions.  Results 
of  the  research  on  these  topics  have  been  reported  to  AFOSR 
in  technical  reports,  and  the  technical  reports  are  included 
as  the  main  body  of  this  final  report. 

Rayleigh  Wave  Dispersion 

Our  studies  using  long  period  data  recorded  in  northeast 
Texas  have  shown  that  a wave  guide  for  Rayleigh  waves  exists 
from  China  and  the  southern  border  of  the  USSR  across  the 
Arctic  to  Texas.  We  determined  group  velocity  dispersion 
curves  for  a number  of  Rayleigh  waves  which  traveled  through 
this  polar  wave  guide  and  found  that  at  periods  greater  than 
30  seconds  the  dispersion  curve  does  not  vary  with  epicentral 
location  within  the  wave  guide. 

Because  of  the  uniformity  of  group  velocity  dispersion 

within  the  polar  wave  guide  and  the  broad-band  nature  of  the 

Rayleigh  waves,  "matched"  or  chirp-filter  techniques  were 

found  to  be  particularly  effective  on  data  recorded  at  our 
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McKinney,  Texas,  site.  A technical  report  to  AFOSR 
entitled  "Continental  Rayleigh  Wave  Dispersion  and  the 
Estimation  of  Chirp-Filter  Detectors"  by  William  Tucker, 

John  McDonald,  and  Eugene  Herrin  is  attached  and  describes 
these  subjects  in  detail. 

Array  Processing  Techniques 
We  have  generated  the  principal  components  of  a major 
soft— ware  system  for  continuous  general  analysis  of  array 
data.  The  system,  called  FKSCAN,  transforms  successive 
blocks  of  array  data  to  the  f raquency-wave  number  domain 
and  explores  that  space  for  correlations.  The  powerful 
advantage  of  transforming  array  data  from  time  and  physical 
space  into  the  frequency-wavenumber  domain  is  that,  quite 
generally,  correlations  are  automatically  separated  and 
sorted  according  to  frequency,  velocity,  and  azimuth 
permitting  detection  and  description  of  signals  not 
readily  discerned  in  time  and  physical  space.  Further, 
spatial  filtering  is  markedly  facilitated  in  transform  space 
and,  as  one  consequence,  large  overriding  signals  may,  in 
effect,  be  "turned  off"  after  the  fact  to  permit  the 
detection  of  much  smaller  simultaneous  arrivals.  FKSCAN, 
after  transforming  and  exploring  each  array  data  block, 
filters  and  removes  the  principal  signals  and  correlations 
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thus  detected  and  searches  the  transform  space  once  again 
for  any  small,  hidden  signals. 

FKSCAN  then  outputs  a bulletin  for  each  data  block  citing 
detections  and  printing  out  the  frequency  power  spectral 
estimates  of  those  signals.  For  each  such  signal  there 
are  also  output  spectra  of  phase-velocity,  back  azimuth, 
and  F-statistic  as  functions  of  frequency.  (The  F-statistic 
is  a measure  of  the  likelihood  that  detected  correlations 
are  genuine;  i.e.,  not  due  to  chance  combinations  of  noise). 

An  auxiliary  program  called  FKPLOT  was  written  which 
outputs  contoured  printer  plots  of  cross-sections  of  frequency 
wavenumber  spectra  cut  normal  to  the  frequency  axis. 


Bias  in  Network  Estimation  of  Body  Wave  Magnitudes 
A method  for  correcting  the  bias  in  body  wave  magnitude 
(mp,)  estimates  is  described  in  an  enclosed  technical  report 
to  the  AFOSR  entitled  "On  the  Estimation  of  Body  Wave  Magnitudes" 
by  Eugene  Herrin  and  William  Tucker.  This  report  shows  that 
observations  from  any  seismolog.ical  network  lead  to  over 
estimation  of  the  magnitude  of  seismic  events  which  are 
near  the  detection  threshold  of  that  network.  Methods  are 
presented  for  calculating  this  magnitude  bias.  Consideration 
must  be  given  to  this  effect  in  comparing  networks  with 
significantly  different  thresholds,  in  comparing  theoretical 
and  empirical  estimates  of  network  capability,  and  in 
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determining  the  source  energy  of  small  seismic  events. 


Effect  on  the  Rayleigh  Spectrum  of  Lateral 
Heterogeneity  in  Earthquake  Source  Areas 

Studies  of  Rayleigh  wave  spectra  obtained  from  a high-gain, 
long  period  seismograph  at  Grand  Saline,  Texas,  show  that 

earthquakes  occurring  in  some  Pacific  island  arc  regions, 
such  as  New  Hebrides,  New  Britain,  and  the  Solomons,  a shift 
of  spectral  energy  tc  longer  periods  with  increasing  source 
depth  is  observed.  This  is  consistent  with  theoretical 
predictions  for  simple  sources  in  laterally  homogeneous  layered 
media.  However,  for  earthquakes  originating  in  certain  other 
island  arcs,  such  as  Tonga— Kermedec , the  Phillipines,  and  the 
Marianas,  no  corresponding  spectral  shift  was  observed.  In  the 
first  group  of  island  arcs,  the  lithospheric  plate  dips  more 
or  less  toward  the  recording  sice  at  Grand  Saline;  in  the 
second  group,  the  dip  is  away  from  Grand  Saline.  It  is  suggested 
that  lateral  heterogeneity  in  the  earthquake  source  regions  is 
at  least  partially  responsible  for  the  spectral  differences 
and  that  in  particular  the  downward  bent  lithosphere  acts  as 
a wave  guide  at  some  azimuths  for  some  Rayleigh  wave  lengths. 

An  enclosed  technical  report  to  AFOSR  entitled  "A  Model  Study 
of  the  Effect  on  the  Rayleigh  Spect  urn  of  Lateral  Heterogeneity 
in  Earthquake  Source  Regions"  by  Tom  Goforth  describes  a 
verification  of  this  hypothesis  using  an  analog  scale  model. 
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The  study  indicates  that  for  intermediate  and  deeo- focus 
earthquakes  in  subduction  zones,  higher  values  of  Ms  will 
be  observed  at  azimuths  opposite  to  the  direction  of  dip 
of  the  lithospheric  plate  than  would  occur  if  the  same 
earthquake  had  occurred  in  laterally  homogeneous  layered 
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Rayleigh  wave  data  were  recorded  by  a single  vertical 

high-gain  seismograph  located  in  north  Texas.  It  was  found 

that  Rayleigh  waves  propagating  in  a continental  path  over 

the  pole  to  this  station  exhibited  very  stable  dispersion 

characteristics , particularly  in  the  period  range  35-75  secs: 

at  shorter  periods  an  average  dispersion  curve  could  be 

estimated.  A composite  dispersion  curve  was  used  to  develop 

"chirp"  (or  matched)  filters,  the  form  of  which  depended  on 

the  distance  between  the  seismograph  and  the  event.  Such 

filters  are  shown  to  be  very  efficient  in  improving  the 

signal-to-noise  ratio  of  Rayleigh  waves  emanating  from  events 

in  the  Sino-Soviet  region  and  can  be  used  as  a means  of 

separating  "mixed"  Rayleigh  waves . For  the  path  over  the 

pole  (or  "polar  waveguide")  the  detection  threshold  is  estimated 

to  be  at  M =3.5. 
s 


VII 


INTRODUCTION 


Ewing  and  Press  (1952)  described  the  classical  "peak- 
and- trough"  method  of  determining  a fundamental  dispersion 
curve  from  a single  surface  wave  coda.  Sato  (1955)  applied 
Fourier  analysis  to  the  calculation  of  dispersion  curves. 

These  techniques  have  been  improved  upon  (Block  and  Hales, 
(1968),  Landisman,  et  al  (1969)  and  Dziewonski,  et  al  (1969)), 
and  an  overall  processing  system  was  devised  combining  Fourier 
analysis  and  time-varying  filters,  enabling  the  calculation  of, 
not  only  the  fundamental  phase  and  group  velocity  curves,  but 
also  the  curves  for  the  higher  modes.  In  particular,  once  a 
group  velocity  curve  has  been  obtained  it  is  possible  to 
design  a time  varying  filter  thus  extracts  a single  dispersed 
mode  and  effects  a significant  improvement  in  the  signal-to- 
noise  ratio.  The  filtered  output  is  then  an  estimate  of  the 
true  signal  Dziewonski,  et  al  (1968)  summarizes  these 
methods  and  it  is  evident  that  they  extract  the  maximum  amount 
of  information  from  the  surface  wave  signal.  Nevertheless, 
the  methods  are  both  time  consuming  and  computationally 
difficult.  But  if  we  restrict  the  investigation  to  the 
fundamental  mode  of  surface  wave  propagation,  it  becomes 


possible  to  develop  a fast,  simple  digital  technique  for 
the  analysis. 

In  this  paper  we  describe  a scheme  which  estimates  the 
fundamental  group  velocity  dispersion  curve  from  a seismo- 
gram, and  then  constructs  a "matched"  filter  based  on  the 
estimated  curve.  This  matched,  or  chirp,  filter  can  then 
be  used  to  operate  on  signals  with  low  signal-to-noise  ratios 


THE  SEISMOGRAMS 


During  the  neriod  March  through  August  1970  a three- 
component  seismograph  system  was  operated  at  a depth  of 
180  m in  the  Morton  Salt  Company  mine  at  Grand  Saline,  Texas 
(Sorrells  et  al.#  1971).  The  seismometers  were  sealed  in 
such  a manner  that  the  masses  were  not  subject  to  buoyancy 
effects  from  atmospheric  pressure  changes.  The  systems  were 
of  the  "advanced"  long  period  type  which  are  operated  with 
a peak  magnification  at  a period  of  50  secs.  Sorrells  et 
al.  (19/1)  have  shown  that  if  such  a system  is  operated  at 
a site  which  is  remote  from  seismic  noise  generated  by  at- 
mospheric processes,  the  spectra  of  the  long  period  noise 
exhibits  a minimum  in  the  period  range  20-60  sec.  Careful 
installation  of  these  seismographs  at  a depth  of  about  180  m 
enabled  the  vertical  seismograph,  the  only  one  used  in  the 
present  study,  to  be  operated  at  a magnification  of  200,000. 
.The  data  were  recorded  on  a digital  acquisition  system  de- 
scribed by  Herrin  and  McDonald  (1971)  after  sampling  at  a 
rate  of  1 per  second. 

The  recorded  seismograms  were  associated  with  source 
parameters  and  body  wave  magnitudes  provided  in  Preliminary 
Determination  of  Epicenters  (PDE)  issued  by  NOAA. 
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DISPERSION  CJRVE  ESTIMATION 


h computer  program  was  written  which  enabled  analog 
representations  of  digital  time  series  to  be  displayed  on 
an  oscilloscope.  Selected  seismograms  were  then  plotted, 
using  a Cal-comp  plotter.  It  was  decided  to  limit  the  ob- 
servations to  Rayleigh  waves  recorded  at  teleseismic  distance 
which  had  traversed  a predominantly  continental  path;  thus 
restricting  the  range  of  azimuths  to  about  15°  east  and  west 
of  the  pole,  and  the  range  of  epicentral  distances  to  55°-130°. 

Using  the  plot  of  the  seimogram  an  estimation  was  made 
of  the  start  time  (Ts)  and  the  length,  in  seconds,  (window) 
of  the  Rayleigh  signal.  A computer  program  was  written  which 
searched  this  window,  starting  at  T , for  zero  crossings,  in 
a manner  similar  to  the  "peak  and  trough"  technique  of  Ewing 
and  Press,  (1952).  The  time  of  the  "actual"  zero  crossing 
was  calculated  by  triangulation  between  the  last  point  immedi- 
ately prior  to  a zero  crossing,  and  the  first  point  immediately 
after  it.  The  initial  zero  crossing  time,  (T  ) was  stored,  and 
at  the  next  zero  crossing,  (T-^ ) the  first  period,  (P^)  was 
computed.  T and  P,  were  then  stored  and  the  seach  continued 
until  the  window  was  exhausted.  For  i = 1,  2,  . . . , N,  the 

group  velocitcs  and  the  periods  were  computed  from 

-2- 
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y = a 

u — ~ 

T'r, 

where  N = the  number  of  zero  crossings  in  the  window, 

A = the  epicentral  distance  from  event  to  station  in 
km,  and 

TT\  = the  observed  travel  time  from  source  to  the  ith 
zero  crossing. 

Since  the  actual  zero  crossing  occurred  within  a one  second 

interval  any  variation  must  be  less  than  one  second.  As  a 

conservative  approximation,  we  assumed  the  error  to  be  uniform; 

2 

thus  the  variance  was  1/12  sec  . 

Now  TT , = T . - 0 
l i 

where  0 is  the  origin  time  of  the  event.  From  Herrin  et  al. 
(1968)  s reasonable  origin  time  variance  is  0.81  to  1.00  sec2. 
Since  lb  and  0 are  independent  observations  their  variances 
may  be  added.  We  thus  take  1.00  sec2  as  a reasonable  variance 
for  each  Tib  . From  Tucker  et  al.  (1968)  and  Herrin  et  al . (1968) 

a reasonable  variance  in  A is  900  km2.  We  then  have 

V.  =.  ^ 

*•  1 > 

T71+ 
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where 


= the  true  distance  (in  km) 


TT , = the  true  travel  time, 

1 

= the  error  in  epicentral  distance  ( (J" ^ = 
= the  error  in  travel  time  ( = 1 sec)  . 

Assuming  an  event  with  = 6000  km  and  a maximum  group 
velocity  of  4.00  km/sec  (Oliver,  1962)  then 


w _ (6000+  ©,)/ 

1 /tisoo+fiz.) 

Since  the  variance  of  e is  small  relative  to  1500  we  can 


ignore  e^  and  write 


Vt  - 4.0  i-  / I5oo 


= 4-.0  S*.  , 


where 


O • 0004-  sec^, 


For  a lower  period  cutoff  of  20  sec  (or  a velocity  of  3.00 


km/sec)  we  get 

VN  =•  (6ooo+  «n)/2ooo 
= 3 o + /2000 

= 3-0  + , 

where 


CT2*^  0’00022S  sec.**. 


Thus  as  a conservative  approximation  we  can  ass 


ume 


= True  velocity  i + £ 


V 


r 2 

where  <JT  = 0.0004  sec  . 

S’ 

Let  us  now  consider  the  variation  in  the  P . 

i 


We  have 


P.  = 2 (T.  -T.  'i 


v.  - I 


whe  re 


O'*  = 0-0833  sec1 
T 


Since  the  minimum  period  is  assumed  to  be  20  sec  the  minimum 
time  difference  is  of  the  order  of  10  sec  and  we  may  assume 
that  the  errors  in  T_^  and  T are  independent.  Then  the 
variance  of  a is 

C1  =4  (O.OS33  r O-OSii)  - O.GGC4  sec1 

Thus  the  variance  in  a velocity  observation  is  three  orders 
of  magnitude  less  than  that  of  a period  observation.  We 
therefore  ignored  the  error  in  the  velocity  computations  and 
assumed  that  the  velocities  were  known  without  error.  Now 
assuming  that  all  the  error  lies  in  the  period  computations, 
it  is  possible  to  consider  a group  velocity  curve  as  a 


regression  of  period  on  group  velocity. 

Figure  2 shows  the  data  points  obtained  from  dispersed 


Rayleigh  waves  for  three  events  located  in  the  Arctic  (Table  1) 
at  similar  epicentral  distances  and  similar  azimuths.  Data 
from  an  event  in  Sinkiang  Province  along  a similar  azimuth, 
but  at  a greater  distance,  (Table  1)  are  shown  in  Finure  3. 

A comparison  of  the  data  in  Figures  2 and  3 in  the  velocity 
range  3.50  to  3.90  km/sec  (periods  approximately  35  and  75  sec 
respectively)  indicates  that  there  is  no  statistical  difference 
between  the  true  dispersion  curve  for  each  event.  The  observed 
scatter  could  result  from  statistical  fluctuations  in  the 
estimates  or  real  differences  in  the  velocity-period  relation- 
ship alang  the  propagation  path.  Thus  it  appears  that  in  the 
velocity  range  3.50  to  3.90  km/sec  there  is  a unigue  disper- 
sion curve  independent  of  path. 

However,  in  the  velocity  range  2.90  - 3.5  km/sec  (periods 
of  20-35  sec),  there  are  indications  that  real  differences  may 
exist  between  the  true  dispersion  curves  in  Figures  2 and  3. 

In  this  velocity  range  we  estimated,  for  the  purposes  of  this 
study,  an  average  dispersion  curve. 

Using  group  velocity  as  the  independent  variable  the 
data  for  each  event  were  grouped  into  cells  of  width  0.5  km/sec. 
Within  each  cell  the  data  were  ordered  according  to  period  value 
and  the  sample  median  was  computed.  The  cell  midpoint 
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velocities  and  median  period:  wore  stored  as  data  pairs. 

The  sample  median  was  chosen  as  a measure  of  central  tendency 
because  it  is  more  efficient  than  the  sample  mean  in  the 
presence  of  outliers  (Dixon,  1953) 

In  the  velocity  range  3.50  to  4.00  km/sec  the  sample 
median  periods  are  observations  on  a single  curve,*  therefore, 
adjacent  points  are  related  and  a smoothing  operator  should 
improve  the  results.  Furthermore,  in  the  velocity  range 
2.90  to  3.50  km/sec  only  a representative  smooth  curve  was 
desired,  and  a three  point  moving  average  operator  was 
applied  to  the  sample  period  values.  The  operator  is  given  by 

Period  (1)  = X(1 ) 

Period  (Num)  = 1/2  • X(Num-l)  + 1/2  • X(Num) 

Period  (i)  = 1/4  • X(i-l)  + 1/2  • X(i) 

+1/4  * X ( i+1 ) 


where  1 

1 

Num 


X ( i ) 
Period  (i) 


2 , . . . , Num-1 

the  index  of  the  cell  with  the  largest 
velocity 

the  index  of  the  cell  with  the  smallest 
velocity 

the  sample  median  period  of  cell  i 
the  smoothed  period  of  cell  i. 
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The  resulting  dispersion  curves  for  the  four  events 
(Figures  2 and  3)  are  shown  in  Figure  4.  The  smoothed  curves 
show  the  same  results  as  the  original  zero  crossing  data. 

That  is,  the  curves  are  not  statistically  different  in  the 
velocity  range  3.50  to  3.90  km/sec,  and,  in  the  range  2.90 
to  3.50  km/sec,  the  smoothed  curves  indicate  that  the  actual 
paths  do  not  produce  markedly  different  dispersions  curves 
at  these  short  periods. 

In  view  of  the  similarity  of  the  curves  we  combined 
the  data  from  the  four  events  and  estimated  a composite 
group  velocity  dispersion  curve  shown  in  Figure  5.  Also  shown 
in  this  figure  is  the  continental  dispersion  curve  compiled 
by  Oliver  (1962) ; the  long  dashed  lines  indicate  bounds  of 
the  expected  scatter  under  normal  observational  conditions. 

The  composite  curve  can  be  seen  to  fall  within  the  error 
bounds.  Furthermore,  the  composite  curve  is  in  good  agree- 
ment with  the  results  of  Ewing  and  Press  (1956)  who  reported 
a group  velocity  of  about  3.9  km/sec  at  a period  of  75  sec. 

The  smoothed  data  points  are  given  in  tabular  form  in  Table  2. 


CHIRP  FILTER  CONSTRUCTION 


We  have  shown  that,  for  our  seismograph  station  in 
north  Texas,  we  can  assume  a single  group  velocity  dispersion 
curve  for  the  continental  path  over  the  pole.  We  can  then 
design  "matched",  or  chirp,  filters  to  aid  in  the  detection 
of  surface  waves  from  small  earthquakes  traversing  the  same 
pa  th . 

A program  was  written  wnich  produced  a sine  wave  that 
decreased  in  period  with  time,  the  dispersion  being  a function 
of  epicentral  distance.  Such  waveforms  have  been  given  the 
name  "chirp"  in  radar  technology.  In  this  case  the  informa- 
tion used  to  distance  - function  were  from  Table  2.  The 
following  assumptions  were  made: 

(i)  the  first  arriving  energy  is  at  the  maximum  period 
observed,  and 

(ii)  between  the  data  points  of  Table  2 the  dispersion 
curve  is  linear. 

With  these  assumptions  we  obtained,  by  simple  algebraic 
manipulations , 

tv  % To=TTu-  Vz; 

v l 


and  for  the  ith_  linear  segment 

_ 2tt  (tfT^V-rrv.,-  1/TTi.) 


-w.  ft)  = 


[P.7  P-  (Pi-,  /it, - O]  ; 


where 


(T.^-T.)  ^ t <(T,-T6)  , 


~yr 

1'  ' = frequency  in  radians  for  the  t th  time, 

= the  i th.  period  in  sec, 

Tit  = the  i th  travel  time  in  sec, 

i - 1 , 2 , . . . , Num. 

The  dispersed  sine  wave  computations  were  given  by 
AMP  (k)  = sin  k t /°u  1 > 

where  k = 1 , 2,  . . . , 

= the  number  of  data  points  in  the  i th  segm  nt 
pi  = phase  shift  of  the  i th  segment 
AMP(k)  = the  dispersed  sine  wave  amplitude. 

For  k = 1 the  segment  was  extrapolated  to  include  time  points 
back  to  Tq.  For  all  other  values  of  k only  time  points  be- 
tween the  TT^  were  included  in  a given  segment. 

The  p . were  as  follows 

Pi 3 0 


Pi 


P 


. - "l/b,  . 

t/  ^ 


( Pi- I 
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is  the  largest  t value  in 


for  1=2,  3,  Num,  where  t 

max 

the  i th  segment.  Also 

t=  k+  TT. 

n-~| 

so  that 


t 

W \i>* 


N^TT..,. 


In  each  linear  segment  the  frequency  was  a function  of 
time  given  by  W^(t).  The  origin  of  the  sine  wave  was  taken 
as  zero,  which  gave  the  time  shift  between  t and  k.  The 
phase  shift  p ^ ensured  a smooth  transition  from  segment  to 
segment . 

In  Figure  6 we  show  the  dispersed  sine  wave  for  one 
of  the  events  used  in  calculating  the  composite  dispersion 
curve  (Figure  5).  These  sinusoids  indicate  tire  effects  of 
differences  in  epicentre]-  distance;  for  the  S ink ia n g Province 
event  (A  = 107.6°)  the  first  zero  crossing  occurs  at  about 
72  sec  and  for  the  other  three  ( ^ = approximately  60°)  the 
first  zero  crossing  occurs  it  about  68  sec.  Similar  period 
differences  are  evident  throughout  the  curves. 

Convolution  of  these  dispersed  sine  waves,  or  chirp 
filters,  w i th  their  respective  signals  should  then  produce 
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psuedo -autocorrelation  functions.  The  results  are  not  true 
autocorrelations  because  the  sine  waves  have  not  been 


amplitude  modulated;  nor  is  the  "true"  dispersion  given 
exactly  by  the  composite  dispersion  curve. 
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APPLICATION  OF  THE  CI1IRP  FILTER 


The  chirp  filter  constructed  from  the  dispersion  curve 
in  Figure  5 was  applied  to  four  different  sets  of  seismograms. 
First,  it  was  convolved  with  the  original  seismograms  used  to 
construct  the  dispersion  curve;  second,  with  other  events 
remote  from  the  first  four;  third,  as  detector  of  surface 
waves  with  a low  signal -to-noise  ratio;  and  fourth,  as  a 
means  of  separating  "mixed"  Rayleigh  waves . 

In  order  to  remove  very  long  period  energy  each  seismo- 
gram was  filtered  prior  to  being  convolved  with  the  chirp 
filter.  The  high-pass  digital  filter  used  had  a corner  at 
50  sec  and  fell  at  36  db/octave  to  a minimum  at  a period  of 
100  sec. 

Figure  7 shows  the  application  of  the  chirp  filter  to 
the  time  series  for  each  of  the  events  used  in  calculating  the 
composite  dispersion  curve.  In  each  example  the  lower  trace 
shows  the  filter  output,  correctly  aligned  with  the  input 
time  series.  The  auto-correlation-like  nature  of  the  output 
should  be  noted;  the  maximum  peak  of  which  should  occur  at  the 
beginning  of  the  Rayleigh  wave  energy.  In  each  of  the  figures 
the  data  were  plotted  so  that  the  largest  excursion  was  full- 
scale;  therefore,  the  improvement  in  signal-to-noise  ratio 
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can  ho  estimated  by  examining  the  reduction  in  noise. 

Three  other  events  were  selected  (table  3)  with  similar 
continental  polar  surface  wave  paths  to  Texas;  however, 
their  epicenters  were  remote  from  thoa?  given  in  table  1. 

The  results  of  convolving  these  seismograms  with  the  same 
chirp  filter  used  previously  are  shown  in  figure  8.  The 
outputs  clearly  show  that  the  filter  is  effective  for  widely 
separated  sources. 

Figure  9 shows  the  results  of  convolving  the  chirp 
filter  with  seismograms  in  which  the  surface  waves  are 
barely  discernible.  This  figure  shows  a seismogram  from  a 
small  event  (m  - 4.7)  at  a depth  of  110  km  in  the  Hindu- 

9 

Kush;  convolution  clearly  emphasized  the  start  of  the  surface 
waves . 

Mixed  Rayleigh  waves  resulting  either  from  multiple 
sources  close  together  or  from  reflections,  or  from  multi- 
pathing  are  very  difficult  to  separate,  even  when  data  from 
large  arrays  are  avialable.  However,  if  the  arrivals  are 
separated  by  times  greater  than  about  50  sec,  the  chirp 
filter  technique  can  be  useful  in  separating  the  arrivals. 
Figure  10(a)  (lower  trace)  shows  the  filtered  output  for  a 
small  event  from  the  region  of  the  USSR-Mongol ia  border; 
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two  and,  possibly,  three  Rayleigh  wave  arrivals  can  be  seen 
on  the  filtered  trace.  The  first  and  largest  of  these 
arrivals  occurs  precisely  at  the  expected  arrival  time  of 
the  fundamental  Rayleigh  wave  for  this  event.  A l airly 
large  earthquake  {fr\  = 5.9)  sparked  an  earthquake  swarm  from 
this  region  of  the  USSR-Mongolia  border  in  May  1970. 

Another  earthquake  from  the  swarm  is  shown  in  figure  10(b); 
again  convolution  detected  two  arrivals.  All  available 
seismograms  of  earthquakes  in  the  swarm  were  processed  but 
only  three  showed  multiple  detections  of  surface  waves. 


DISCUSSION 


We  have  described  above  a new  approach  to  the  "peak 
trough  method  of  surface  wave  analysis  which  has  been 
used  to  determine  group  velocity  dispersion  curves.  The 
machine  measurements  proved  to  be  at  least  as  good  as  those 
of  an  experienced  analyst,  and  could  be  produced  much  faster. 

A further  advan’age  was  gained  in  the  present  analysis 
from  the  special  properties  of  the  seismograph  system. 

Many  previous  surface  wave  studies  have  used  data  from  the 
World-Wide  Standard  Seismograph  Network  recordings  in  which 
the  peak  period  of  the  calibration  response  is  at  20  sec. 

In  our  systems,  the  peak  response  proved  to  lie  within  the 
"stable "portion  of  the  group  velocity  dispersion  curve. 

This  type  of  analysis,  while  being  quick  and  computa- 
tionally easy,  nevertheless  has  its  problems.  Any  non-least- 
t ime  arrivals  will  impair  the  estimation  of  the  dispersion 
curve?  the  technique  does  not  give  as  much  information  as 
that  of  Dziewonski  et  al . (1969),  which  places  a restriction 
or  the  use  of  the  method. 

Nevertheless,  by  carefully  selecting  events  with  a 
large  signal-to-noise  ratio,  and  with  an  absence  of  either 
higher  modes  or  very  long  period  surface  waves,  it  has  been 
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shown  that  a single  fundamental  group  velocity  dispersion 
curve  can  be  calculated  for  a continental  path  over  the 
pole  to  the  central  United  States.  This  claim  can  be  further 
substantiated  by  the  fact  that  a chirp  filter,  constructed 
from  this  fundamental  dispersion  curve,  proved  effective 
for  events  from  a region  extending  from  southern  USSR  and 
northern  China  to  the  Arctic  Ocean. 

In  addition  to  the  events  mentioned  in  the  previous 
sections  of  this  paper,  seismograms  from  the  other  events 
listed  in  table  3 were  processed.  In  some  cases  the  auto- 
correlations displayed  low  s ignal-to-noise  ratios,  but  in 
every  search  for  an  event  reported  by  NOAA,  an  event  was 
detected;  unreported  events  were  also  detected.  Thus  further 
study  is  required  to  determine  operational  detection  thres- 
holds and  false  alarm  rates. 
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CONCLUSIONS 


In  view  of  the  results  reported  here  it  might  be  appro- 
priate to  suggest  a reason  for  the  stability  of  the  group 
velocity  dispersion  curve,  particularly  at  the  longer  periods, 
on  a continental  path  over  the  pole. 

It  has  been  shown  in  recent  years  (e.g.  Roy  et  al.  1972) 
that  the  western  United  States,  known  structurally  as  the  Basin 
and  Range  Province,  is  characterized  by  high  temperatures  in 
the  lower  crust  and  upper  mantle.  This  region,  which  is  known 
to  extend  up  the  western  part  of  Canada,  is  also  a poor  trans- 
mitter of  long  period  surface  waves.  A particularly  good 
example  of  this  phenomenum  can  be  seen  in  figure  11.  The 
Greenland  Sea  event  {table  1),  as  recorded  at  Grand  Saline, 
is  compared  to  the  same  event  recorded  at  Queen  Creek,  Arizona. 
This  station  (Fix  and  Sherwin,  1970)  is  located  in  the  heart 
of  the  Basin  and  Range  Province.  The  gains  of  the  seismographs 
at  the  two  stations  were  equalized  at  50  sec,  but  it  will  be 
seen  that  the  longer  period  surface  waves  art  virtually  miss- 
ing at  the  Queen  Creek  site. 

To  the  east  the  continental  United  States  are  characterized 
by  the  Atlantic  and  Gulf  coastal  plains.  Furthermore,  the 
central  plains  are  known  to  have  relatively  low  temperatures 


in  the  lower  crust  and  upper  mantle  (Roy,  et  al.,  1968, 

Combs  and  Simmons,  1973)  which  provides  an  exceptionally  good 


path  for  Rayleigh  waves  with  periods  from  20-75  sec.  Our 
detection  station  in  north  Texas  is,  essentially,  at  the 
southern  end  of  a surface  wave  guide  which  passes  over  the 
north  pole  (Figure  12)  . 

The  analog  plots  of  the  seismograms  used  in  the  study 
were  analysed,  after  they  had  been  high-pass  filtered,  for 
the  presence  of  detectable  Rayleigh  waves.  Using  the  formula 

M = log  + 1.66  log  A - 0-18,  which 

can  be  easily  derived  (Gutenberg,  1945),  values  of  Ms  were 
found  where  A is  the  penk-to-peak  amplitude,  in  mp,  of  the 
surface  wave  at  a period  of  T=20  sec,  and  A is  the  distance 
from  epicenter  to  detector,  in  degrees.  The  results  of  these 
computations  are  given  in  tables  1 and  3,  where  the  values  of 
Ms  are  compared  to  the  values  of  m^  determined  by  NOAA. 

Using  the  data  presented  in  this  study  an  interim  detec- 
tion threshold  has  been  determined  for  our  Grand  Saline  record- 
ing station.  For  events  travelling  within  the  polar  waveguide 
this  threshold  is  at  Ms  = 3 *5,  for  a single  vertical  seismometer, 
which  is  equivalent  to  a threshold  of  m^  = 4*5  for  events  at 
depths  of  less  than  50  km. 
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FIGURE  CAPTIONS 


1.  Response  of  the  seismograph  system. 

2.  Rayleigh  wave  dispersion  information  recorded  in  north 

Texas  for  three  events  (see  table  i),  (a)  North  of 

Svalbard,  (b)  Greenland  Sea,  (c)  Norwegian  Sea. 

3.  Rayleigh  wave  dispersion  information  recorded  in  north 
Texas  for  an  event  in  Sinkiang  Province  (see  table  1) . 

4.  Rayleigh  wave  dispersion  curves  for  the  four  events  in 
table  1 after  smoothing  (a)  North  of  Svalbard  (b)  Green- 
land Sea  (c)  Norwegian  Sea  (d)  Sinkiang  Province. 

5.  Composite  Rayleigh  wave  dispersion  curve  for  the  four 
events  in  table  1 in  conjunction  with  continental  disper- 
sion curve  due  to  Oliver  (1962) . 

6.  Typical  dispersed  sine  wave,  or  chirp  filter,  derived 
from  the  dispersion  curve  shown  in  figure  4 (b) . 

7.  Convolution  of  the  chirp  filter  derived  from  the  curve 

in  figure  5 with  each  of  the  four  events  given  in  table  1, 
(a)  North  of  Svalbard,  (b)  Greenland  Sea  (c)  Norwegian 
Sea  (d)  Sinkiang  Province. 

8.  Convolution  of  the  chirp  filter  derived  from  the  curve 
in  figure  5 with  three  of  the  events  given  in  table  3: 
events  from  very  different  locations. 

9.  Convolution  of  the  chirp  filter  derived  from  the  curve  in 
figure  5 with  an  event  from  the  Hindu-Kush  (table  3)  in 
which  surface  waves  are  barely  discernible. 

10.  Convolution  of  filter  and  surface  waves  showing  two,  or 

possibly,  three  surface  wave  arrivals  (a)  USSR  Mongolian 
border  on  17  May  1970,  - 4.5,  (b)  USSR  - Mongolian 

border  on  23  May  1970,  = 4.5.  (See  table  3 for  details). 


23< 


11.  Comparison  of  the  surface  waves  from  the  Greenland 
Sea  event  (table  1)  recorded  at  (a)  Queen  Creek, 
Arizona,  (b)  Grand  Saline,  Texas. 


12.  A polar  plot  of  the  world  centered  on  Grand  Saline. 
The  locations  of  the  epicenters  of  the  events  used 
in  this  study  are  shown;  the  alphanumeric  designation 
refers  to  tables  1 and  3. 


TABLE  CAPTIONS 


1.  Data  for  the  four  events  used  to  calculate  the 
composite  dispersion  curve  shown  in  figure  5. 

2.  Coordinates  of  the  points  used  to  plot  the  curve 
in  figure  5. 

3.  Data  for  the  other  events  used  in  this  study. 


onmiN  • o 1:0  r.  r q io, 
twin  i rr  suuiJiino 


•O  tij  A)  80  'JO  HD 


Figure  2 


□mem  ■ o 139  ? ? ii. 
CH.Lii  nm  srn 


•jcT  to  vu  iu  ao  iuq 


11  lll«) 


4.25 


mmm 


vi:\(Aix  > 

tyco'jo.  _>iu  .y.x> 


ORIGIN  • 0 210  5 
SlSKIftNC  PROUIfCE 


*•  • rii 


3- 1*3- l*  31- 


CEL-St.9  **M1. 


TABLE 


to 

U 

P 

ro 

u 


w 

Q 

Cu 


C 

< 

O 

£ 


E 

o 

p 

tn 


n 

z 

Z 

0 

cr> 

O 

■P 

M-t 

< — t 

rH 

ro 

o 

rH 

- 

- 

P 

T3 

CP 

CM 

U 

c 

m 

rH 

ro 

(0 

O 

o 

u 

s-i 

(N 

in 

0 

m 

tr> 

•P 

C 

c 

•H 

0 

10 

to 

u 

ro 

P 

•H 

X 

a 

tu 

— . 

0 

Eh 

in 

■^r 

OH 

QJ 



cp 

O 

C 

* 

rH 

•rH 

m 



to 

rH 

0 

ro 

0 

CP 

-P 

c/j 

rH 

P 

to 

0 

c 

tj 

ro 

rO 

•H 

C 

■p 

c 

T> 

ra 

0 

P 

p 

o 

V 

-P 

p 

O 

O 

O 

jj 

cn 

o 

u 

c 

>1 

ivj 

o 

rH 

E 

E 

*H 

Or 

o 

o 

-p 

n< 

P 

P 

ro 

ro 

OP 

U-f 

-p 

10 

to 

X> 

TJ 

0 

0 

0 

Cp 

-p 

-P 

-P 

c 

o 

ro 

ro 

•H 

c 

rH 

rH 

XJ 

-p 

D 

D 

p 

o 

U 

O 

0 

o 

rH 

»H 

u 

0-1 

ro 

03 

0 

u 

U 

p 

0 

E 

ro 

• 

• 

U) 

CN 

ro 

TABUS  2 


Gr.  ou 


p Vrloci  L\' 

Period 

3.9250 

75.5223 

3.3750 

69.2352 

3.3250 

60.8679 

3.7750 

51.3606 

3.7250 

45.6151 

3.6750 

41 .9743 

3.6250 

39.5603 

3.5750 

37.4096 

3.5250 

35.5422 

3.4750 

34.0502 

3.4250 

32.0733 

3.3750 

30.5458 

3.3250 

29.3486 

3.2750 

27.7666 

3.2250 

26.1123 

3.1750 

24.6726 

3.1250 

23.4360 

3.0750 

21.3606 

3.0250 

19.2884 

2.9750 

18.6214 
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(For  footnote  see  table  1) 
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ABSTRACT 


Observations  from  any  seismological  network  lead  to 
overestimation  of  the  magnitude  of  seismic  events  which 
are  near  the  detection  threshold  of  that  network.  Methods 
are  presented  for  calculating  this  magnitude  bias.  Con- 
sideration must  be  given  to  this  effect  in  comparing  networks 
with  significantly  different  thresholds,  in  comparing 
theoretical  and  empirical  estimates  of  network  capability 
and  in  determining  the  source  energy  of  small  seismic  events. 
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INTRODUCTION 

Body  wave  magnitude  was  first  defined  by  Gutenberg 
(1945  a,b),  and  his  distance  and  focal  depth  normalization 
factors  are  still  in  general  use.  Recently  Veith  and  Clawson 
(1972),  based  on  a large  quantity  of  data  from  explosions  and 
earthquakes,  have  revised  these  normalization  factors.  In 
this  paper  we  assume  that  the  calculations  required  to  obtain 
magnitudes  have  been  made  correctly,  that  the  depth-distance 
corrections  are  exact  and  that  there  are  no  gross  errors  in 
the  observations.  We  are  primarily  concerned  with  problems 
which  arise  in  the  determination  of  magnitudes  for  small  seismic 
events  which  are  near  the  detection  limit  of  the  observing  network. 

Freedman  (1967)  showed  that  the  individual  observations  of 
amplitude  divided  by  period  (A/T)  used  in  the  estimation  of 
magnitude  for  a given  event  tend  to  be  log-normally  distributed; 
thus,  the  values  of  log  A/T  are  approximately  normally  distributed. 
She  concluded  that  the  mean  of  the  magnitudes  determined  by  each 
station  should  be  used  to  characterize  the  source  strength  of  a 
given  event.  Freedman  worked  with  large  events,  however,  and 
was  not  concerned  with  the  problem  of  missing  data.  In  the 
following  discussions  we  define  the  "true"  magnitude  of  an 


event  to  be  the  mean  of  the  individual  magnitudes  as  observed 


by  stations  in  a "perfect  network";  a "perfect  network" 
for  which  the  probability  of  each  station  observing  the 
from  the  event  is  essentally  unity. 
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DISTRIBUTION  OF  OBSERVED  MAGNITUDES 

Mack  (1969)  showed  that  surprisingly  large  variations  in 
amplitude  were  observed  among  the  elements  of  LASA  for  the 
first  arriving  P-waves  from  any  given  event.  He  attributed 
this  variation  to  a multipath  phenomenon.  We  have  investi- 
gated the  variation  in  magnitude  estimates  which  would  arise 
from  the  multipath  effect  proposed  by  Mack. 

We  assumed  an  observed  signal  to  be  S(t)  where 


5 (- 1)  = Si  ( f)  + dS,{  t +TI)  4-  b S»  (t  +T3) 


(l) 


S^t)  is  the  LASA  summation,  properly  phased,  for  a presumed 
explosion  on  April  21,  1966,  in  East  Kazahkstan,  U.S.S.R. 
(figure  1).  S(t)  was  computed  and  plotted  for: 


a 

= 0. 
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b 
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r 3 

= 0, 

0, 

.1, 

0.2,  . 

d plotted 

2 U)  = S,  (i)  4-  c 5,  Th) 


(2) 
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for 


c 


0.2,  0.4,  ...  , 1.0  and 


^4  0.1,  0.2,  ...  , 1 . 8 sec. 

These  combinations,  plus  SL(t),  provided  2596  signals.  The 
peak-to-peak  amplitude  (A)  and  the  period  (T)  were  read  for 
the  maximum  peak-to-peak  excursion  in  the  first  5 sec.  of 
each  signal,  according  to  the  method  used  by  analysts  report- 
ing LRSM,  WMO,  BMO,  CPO,  UB0  and  TFO  data.  Histograms  for 
the  observed  values  of  A/T  and  log  A/T  are  shown  in  figures  2 
and  3.  These  histograms  should  be  compared  with  figures  1 
and  2 in  Freedman  (1967).  The  variation  in  signal  amplitude 
seen  in  this  study  simulates  the  effects  to  be  expected  from 
multipath  phenomena  of  the  kind  suggested  by  Mack  (1969) 

The  variations  of  A/T  tend  to  be  log-normally  distributed: 
multipath  effects  may,  in  part,  explain  the  empirical  results 
reported  by  Freedman  (1967). 

In  the  following  sections  we  have  assumed  that  the 
individual  observations  of  magnitude  for  an  event  are  normally 
and  independently  distributed.  In  our  simulation  study  the 
standard  deviation,  for  2596  observations  of  magnitude,  was 
0.12  magnitude  unit.  Freedman  (1967)  reports  a standard 
deviation  (SD)  for  USCGS  data  of  about  1/2  unit,  and  Mack 
(personal  communication,  1972)  reported  a SD  for  hASA  of  about 


fl 
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1/3  unit.  Veith  and  Clawson  (1972)  give  values  of  SD  of 
0.356  for  explosion  data  and  0.350  for  earthquake  data. 

We  believe  that  the  simulation  studies  give  a low  value  of 
amplitude  variability  because  only  the  effects  of  three  com- 
binations of  a simple  P-wave  signal  were  considered.  Thus 
reasonable  limits  for  the  SD  of  the  magnitude  distribution 
are  0.2  to  0.5  unit;  a value  of  0.35  magnitude  unit  has  been 
used  for  the  standard  deviation  in  all  calculations  reported 
in  this  paper. 
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ESTIMATION  OF  MAGNITUDE  FOR  SMALL  EVENTS 


In  any  given  network  various  stations  will  fail  to  detect 
and  determine  magnitudes  for  smaller  events.  We  speak  of  a 
90%  body  wave  detection  threshold  as  that  magnitude  for  which 
there  is  a probability  of  0.9  that  at  least  n stations  of  the 
network  will  observe  the  initial  P-wave,  but  only  those  stations 
observing  the  signal  can  report  a magnitude. 

Stonewall  model  - Figure  4 (a,  b,  c)  shows  a simple  model  which 
illustrates  the  problem  of  estimating  magnitudes  for  small 
events.  In  this  model  no  station  can  observe  a magnitude  less 
than  the  stonewall  value,  m^  . As  long  as  the  "true"  magnitude 
of  the  event,  m^,  is  large  with  respect  to  m^  (figure  4a),  the 
expected  value  of  the  observed  magnitude  is  the  "true"  value. 
When  approaches  md  (figure  4b)  only  those  values  of  magni- 
tude to  the  right  of  the  s tonew all  can  be  observed  by  the 
network,  thus  the  expected  value  of  the  observed  magnitude, 

E(m),  is  larger  than  the  "true"  value.  The  resulting  bias 

increases  with  decreasing  values  of  m . An  extreme  case  is 

b 

shown  in  figure  4c  where  the  "true"  magnitude  is  less  than 
the  stonewall  value.  Notice  that  the  reported  magnitude; 
i.e.,  the  mean  of  the  observed  magnitudes,  can  never  be  less 
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tAe 

thanA stonewjll  magnitude,  no  natter  how  small  the  "true"  value 
becomes . 

The  s tonewall  model  is  an  obvious  simplification  of  what 
really  happens  in  a network,  but  it  can  be  used  to  illustrate 
the  effect  with  which  we  are  concerned.  The  model  describes 
an  automated  detection  system  using  signals  from  a large  array. 
In  this  case  the  distances  from  source  to  each  element  of  the 
array  are  nearly  the  same  and  noise  levels  should  not  vary 
significantly  among  the  elements.  If  a fixed  threshold  is 
used  in  the  detection  program,  then  the  stonewall  model  is 
an  accurate  description  of  the  process. 

Rar.dom-wall  model  - A more  useful  model  for  the  study  of 
magnitude  bias  is  one  in  which  the  location  of  the  stonewall 
is  a random  variable;  that  is,  we  can  only  speak  of  the  prob- 
ability of  the  wall  having  a particular  value  for  any  given 
observation.  The  variability  in  location  of  the  wall  depends 
on  such  factors  as  differences  in  epicentral  distances,  in 
noise  levels  and  i.n  detection  criteria.  It  is  impractical  to 
determine  a.  priori . the  distribution  function  for  the  random- 
wall  . We  have  carried  out  the  computations  reported  in  the 
next  sections  for  two  distributions,  the  nocmal  and  the  uniform. 
Conclusions  are  essentially  the  same  for  both  cases;  all  of 


} 
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the  following  equations  and  calculations  assume  a normal 
distribution  function  for  the  random-wall . 

Terms  are  defined  as  follows: 

booy  wave  magnitude  - a normal  random  variable 
mean  of  - "true"  magnitude 
standard  deviation  of 

Random-wall  - a normal  random  variable 
mean  of  Md 

standard  deviation  of  M, 

d 

2 2 

= cr  + <r 

b d 

is  the  cumulative  F - distribution  with  k,l 
degrees  of  freedom 

is  the  probability  that  at  least  n stations  from 

n,  N 

a total  of  N stations  detect  an  event  with  "true" 

magnitude  m^. 

A 

m magnitude  observed  by  those  stations  that  detect 
the  event  - a random  variable 
We  define  p as  the  probability  of  detection  of  an  event 
by  an  individual  station;  thus 
p = Pr  (Mb  > Md)  . 

Letting  i = 2n  and  j = 2 (N-n  + 1)  we  have  that  (see  appendix 


*b 

™b 

tf~b 

Md 

m. 


<r 


cr 


k.  1 
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for  derivations) 


P = f / • 
n , N ' '-,J 


and 


E(m  nt  m)  x m 


± -t 

I - -p> 


(3) 


(4) 


The  LHS  in  equation  (4)  is  the  conditional  mean  of  m given 
that  at  least  n of  N stations  in  the  network  detected  the 
event.  The  quantity  E(m  |n,N)  is  defined  as  the  apparent 
magnitude  of  the  event;  the  magnitude  we  would  expect  the 
network  to  report. 
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NUMERICAL  RESULTS 


Computer  algorithms  were  written  using  equations  (3)  and 
(4).  The  first  code  gives  the  detection  probability  as  a 
function  of  "true"  magnitude.  In  equation  (4)  we  see  that 
the  second  term  on  the  RHS  is  the  magnitude  bias.  The  computer 
code  calculates  the  apparent  or  biased  magnitude  as  a function 
of  "true"  magnitude. 

In  order  to  arrive  at  reasonable  estimates  of  the  stan- 
dard deviation  (CT^)  of  M^ , the  random-wall  distribution,  we 
must  know  something  about  the  properties  of  the  network.  If 
we  restrict  epicentral  distances  to  lie  between  30  and  70  degrees 
the  range  of  log  A/T,  from  the  distance  effect  alone  is, 
according  to  Veith  and  Clawson  (1972),  about  0.2  magnitude 
unit.  Noise  at  individual  stations  may  well  change  in  ampli- 
tude by  a factor  of  two  giving  a variability  in  the  detection 

threshold  of  about  0.3  magnitude  unit.  Clearly  when  = o 

d 

we  revert  to  the  stonewall  model . For  the  restricted  network 
case  we  take  = 0.20,  a value  which  appears  to  allow  fcr 

the  expected  variability.  Tables  1 and  2 and  figures  5 and  6 
show  calculations  for  two  networks.  The  first  is  character- 
istic of  the  existing  WWSSN  with  its  90%  detection  threshold 
reported  to  be  about  4.5  magnitude  units.  The  second  is  a 
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An  event  which  appears  to  the  WWSSN  to  have  a magnitude 
of  4.5  has  a "true"  value  of  4.1.  This  same  event  will  appear 
to  the  Improved  Network  to  have  a magnitude  of  4.2.  Both 
estimates  are  biased;  the  relative  bias  is  0.3  unit. 

Figures  5 and  6 show  the  interval  probability  distribu- 
tion (see  SIPRI,  1968,  p.  48)  for  the  two  networks.  The  spread 
in  apparent  magnitude  between  the  10%  and  90%  is  only  about  0.1 
unit,  smaller  than  observed  for  real  networks.  Studies  of  de- 
tection capability  for  networks  composed  of  LRSM  stations, 

WWSS  and  stations  such  as  TFO  (see,  for  example,  Tech.  Rept. 
64-163,  VELA  Project  VT/036,  1964,  Geotech -Teledyne)  show  that 
the  apparent  range  betv een  10%  and  90%  detection  probability 
is  more  like  0.3  or  0.4  magnitude  unit.  Of  course,  in  the 
real  network  studies  there  is  no  restriction  on  epicentral 
distance;  all  detected  events  are  included.  As  predicted  by 
Evernden  (1967)  and  observed  by  Veith  and  Clawson  (1972),  the 
variability  of  signal  levels  is  greatly  increased  when  the 
entire  distance  range  for  direct  P arrival  is  considered. 
Because  we  do  not  know  a.  priori  what  the  epicentral  distances 
will  be  for  our  hypothetical  networks , we  must  model  the 
effect  of  removing  the  distance  restriction  by  increasing  the 
standard  deviation  of  md . We  found  that  with  O-,  about  0.5 
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the  apparent  spread  in  magnitude  from  a 10%  to  90%  detection 
probability  is  about  0.35  unit,  in  good  agreement  with  ob- 
served results.  Table  4 and  5 and  figures  7 and  8 show 
detection  capability  and  magnitude  bias  for  the  two  hypothe- 
tical networks  with  <5“d  = 0.47.  A summary  of  these  results  is 
given  in  Table  6. 

TABLE  6 

Unrestricted  Distance  Model 


Detection 

Apparent 

True 

Probabilitv 

Maqnitude 

Maqnitude 

WWSSN 

90% 

4.50 

4.25 

50% 

4.33 

4.03 

10% 

4.14 

3.77 

Improved 

Network 

90% 

3.89 

3.65 

50% 

3.72 

3.41 

10% 

3.54 

3.17 

From  these  calculations  we  see  that  an  event  which  appears 
to  the  WWSSN  to  have  a magnitude  of  4.5  has  a "true"  value  of 
4.25.  This  same  event  will  appear  to  the  Improved  Network  to 
have  a magnitude  of  4.36.  Again,  both  estimates  are  biased; 
the  relative  bias  is  0.14  unit. 
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DISCUSSION 


Simulation  studies  have  been  used  to  determine  the 
theoretical  detection  thresholds  for  various  networks  (see 
SIPRI,  1968,  pp.  115  - 125).  Maps  produced  from  these  studies 
show  contours  of  detection  threshold  in  terms  of  "true"  mag- 
nitude, as  defined  earlier  in  this  paper.  Empirical  studies 
of  data  obtained  from  one  of  these  networks  would  provide 
detection  thresholds  in  terms  of  apparent  magnitude.  Unless 
magnitude  bias  is  considered,  no  real  network  can  ever  appear 
to  have  as  low  a detection  threshold  as  was  predicted  by  the 
simulation  studies.  The  computer  codes  used  for  these  studies 
can  easily  be  modified  to  give  the  apparent  magnitude  for 
each  event  detected;  all  of  the  pertinent  information  is  a- 
vailable  in  memory.  We  suggest  that  network  studies  of  this 
kind  should,  in  the  future,  provide  threshold  calculations  in 
terms  of  both  apparent  and  "true"  magnitude,  and  that  the  re- 
sults of  empirical  network  studies  should  be  corrected  for 
magnitude  bias . 

The  computations  of  magnitude  bias  for  small  events  pre- 
sented in  the  previous  section  are  conservative  in  that  we 
assumed  an  accurate  measurement  of  the  signal  amplitude  had 
been  made;  additive  noise  was  not  considered.  In  general, 
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the  presence  of  noise  will  increase  the  magnitude  bias  for 
threshold  events.  Results  for  the  two  networks  show  conclusively 
that  magnitude  bias  must  be  taken  into  account  in  the  comparison 
of  performance  by  networks  with  significantly  different  thres- 
holds and  in  the  estimation  of  the  source  energy  of  events  near 
the  threshold  for  the  detecting  network. 
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APPEMDIX 


Derivation  of  the  net:  detection  probability.  We  suppose  that 
the  observed  body  wave  magnitude,  is  a random  variable  with 

Mb  - M (w\  b,  ) 


and  that  the  threshold  magnitude,  M . is  also  a random  variable 

d 


with 


Mol  ~ M (MJ,  ) 


with  and  independent.  Then  the  probability  that  any  given 
station  detects,  p,  is  given  by 

jo  = Pr( Mb>/Vtd)=  Phi'Mk-McJ  >o)  _ 


But 


Mk-Md 


where 


so  that 


Z ■z-  z' 

<r  = -t 


o / Mb  - Md  -(Mb-Yfld) 

sp  = r h(  


> 


md  - Wb 


= Pt(z 
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Krtd  - nib 


C5 


— 1 " 


/v-i 


< 


A'  'i- 


- 


where  Z is  the  standard  normal  variable.  Thus  we  obtain  the 
result  that  the  individual  station  detection  probabilities  are 
all  equal.  It  then  follows  from  the  definition  of  the  binomial 
distribution  that 

Ph  ( N ct  Detect  i on')  = P h f N 


N-  x 


where 


n\ 


r 

N l 


x / " (AJ-x) ! X J 

But  as  is  well  known  the  binomial  cumulative  can  be  evaluated  by 


N - 


d X 


so  that 


and 


d X - 
d 2 

Pm,  n = 


- 2 


(/  + 2) 
r(N+i) 
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/ -p 


^ (i 


z)  dz 


which  gives 
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J = Z-(  N - 

Now  let 

=t  = T * / 0 

so  that 

J z-  -4- 

ol  t J 

and 

£ + J 

2. 


r »-p 


■ ±-i 

Z z 
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«+*)  i; 
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r^) 
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/ - P 
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which  is  the  cumulative  F.  Thus 


where 


p,/N-  Fi,j(4£) 


l = 2 n 

j - Z ( N-  h -h  l) 


P 


= Pr{  ? 


> 


/p7o/  - m b 


We  note  that  the  result  is  valid  for  any  constant  station  detection 
probability,  p,  so  that  other  distributional  assumptions  on  and 
are  possible.  Note  that  if  CPJ  - o then  degenerates  to  the 
point  m^  and 

P 


- P,[z  > 
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giving  the  stonewall  case.  W - the  other  parameters  held  fixed 


the  effect  of  increasing  CS^  is  to  rotate  p about  m^  = m^ . Thus 
f°r  md  p decreased  and  for  m^<  m^  p is  increased.  If  p 

p 

increases  then  increases  even  more  and  if  p decreases  ~~I~p 


/-/=> 

decreases  even  more.  Therefore,  with  the  other  parameters  fixed 

the  final  effect  on  P is  a rotation.  For  m > m F is 

n , N b d i ' j 

decreased  as  csQ  is  increased  and  for  < n»d  Fi  is  increased 
as  CQ  is  increased. 

It  is  easy  to  show  that 

p I , /Mb 


•/  1 


and 


P 


— ^>  o , (X.  5 yyi  t>  — > — <=►<> 


If  m^  is  changed,  say 


and  we  let 


rr)d  - + a 

nib  - Wb  + A 


then  p remains  the  same  ( assuming 'a_"  to  be  constant)  . Thus  with 
the  other  parameters  held  fixed  an  incremental  change  in  m^  pro- 
duces only  a location  shift  in  P 

n,  N 

Derivation  of  the  bias  in  magnitude  estimates.  We  suppose  that  the 
station  detection  probability,  p,  is  constant  and  that  the  magnitude 

A 

estinate,  m,  is  formed  from  only  those  stations  which  detect. 


Thus 


A 

m ~ 


nd 
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where 


n = the  number  of  st  itions  that  detect, 
d 

M^(k)  = the  magnitude  of  the  k-—  station  that  detected. 
Since  only  the  detected  magnitudes  are  employed  any  expectations 
must  be  formed  conditional  on  a detection  occurring  and  the  mean 

A A I A I 

of  m is  E (m  I Net  dection)  = E(m  n,N),  say.  We  have 


A|h,n) 


With  our  previous  assumptions  on  and  we  have  that 
E(Mb(k)  j M^k)  >Md(k))=  E(Mh  j Mb  > M ) t 


Then 


E (m  n,  N) 


K>d  = n 

= E [m  b ) Mb  > Md] 


Prt.N 


Thus  the  mean  of  the  magnitude  estimator  is  a function  of  only  the 

conditional  distribution  of  given  . This  result  also  holds 

for  any  distributional  assumptions  on  and  that  give 

E(M  (k)  M (k)  >M  (k ) ) a constant  over  k and  p constant, 
b I b d 

Let  the  conditional  density  be  denoted  by  f <*l  Hb>  M ) • 

Then 

f (x  I Mb>Md)  * 


-5- 


V J' 


d X 


Ph  ( Mb  - ^ Mb  ^ ^ d ) 


Pr  ( Mb  ?Md) 


~ -i-  FrfMb  <X  aJ  Mb>Md) 
r d x ' 


We  have 
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We  then  have  that 
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The  conditional  mean  is  thus 
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But  the  first  integral  is  just  p flff'  so  that 
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Thus  A4b>Mdj=  ^V1fc>  + 0 <s~ 


upon  completing  the  square  in  the  exponent. 
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where 
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o~ Z-  ° ■+■  o; 
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CS~t, 
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IfcSd  = 0 then  degenerates  to  the  point  m^  and 


E(m  ) n,N)  = KV^  b + ~p£rr 
giving  again  the  stonewall  case. 

Recall  that  as  m^— > °*>  , p — ? I . In  this  case  it  is  easy 
to  show  that 

E(m  | n,N) > Mb,  Mb  > 

In  the  case  that  m^->_00,  p ° and  we  employ  L'Hopital's  rule 


to  show  that 

E (m  | n,N) 
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Thus  if  = O , then 


(m  n,  N) 


— > m 0^3  VY\  b 

d ' 


which  is  intuitively  clear.  If  then 


E (m  i n , N) 


cr  2 / 


0.5  n — 


and  E(m  ) n,N)  is  asympototic  to  . In  this  case  the  effect 
is  a reduction  in  bias  relative  to  the  case  when  = 0 for 

"smalle;"  m^. 

Other  distributional  assumptions.  We  also  studied  the  case  of  a 
uniform  distribution  on  and  a normal  on  M^.  That  is  we  supposed 

Mb  ~ N ( Mb,  oV) 

Md  \J  ( Wd  J VJidih  Wd) 

and  M and  M,  independent.  After  straight  forward  but  tedious 
Vi  0 


computations  one  obtains 
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Also  one  has 
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and 
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It  is  then  easy  to  show  that 
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The  case  in  which  - 


t . CLS  KHl 


<=>o 


— oo  is  not  as  straight  forward 


However,  in  this  case  - w^  is  a lower  bound  on  any  magnitude 


that  is  observed  with  a nonzero  probability.  Thus  in  this  case 


E(Mh\Mh>  Md)  ^ md  - wd 


for  all  m.  . Relative  to  the  normal  distribution  case  the  bias 
should  be  greater  for  smaller  values  of  np  . Test  cases  verify  this 
In  fact  for  comparable  values  of  c2b,  and  Cy  and  w^  the  detection 
probability  curves  were  essentially  equal.  In  every  case  at 


values  of  mfa<  md  the  bias  in  the  uniform  case  was  greater  than 


that  in  the  normal  case. 
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ABSTRACT 


MM 


Theory  has  been  previously  developed  to  determine  the 


spectrum  of  Rayleigh  waves  which  have  been  excited  and  propa- 


gated through  laterally  homogeneous  layered  media.  The  spectral 


shape  is  a complicated  function  of  the  source  and  medium  proper- 


ties, but  in  general  there  is  a shift  of  spectral  energy  in 


relative  sense  toward  the  longer  periods  as  the  source  depth 


is  increased. 


Studies  of  Rayleigh  wave  spectra  obtained  from  a high-gain, 


long  period  seismograph  at  Grand  Saline,  Texas,  show  that  for 


earthquakes  occurring  in  some  Pacific  island  arc  regions,  such  as 


the  New  Hebrides,  New  Britian  and  the  Solomons,  the  relative  spec- 


tral shift  to  longer  periods  with  increasing  source  depth  is  indeed 


observed.  However,  for  earthquakes  originating  in  certain  other 


island  arc  regions,  such  as  Tonga -Kermedec,  The  Phillipines,  and 


the  Marianas,  no  corresponding  spectral  shift  was  observed.  In 


the  first  group  of  island  arcs,  the  lithospheric  oiate  dips  more 


or  less  toward  the  recording  site  at  Grand  Saline;  in  the  second 


group,  the  dip  is  away  from  Grand  Saline.  It  is  suggested  that 


lateral  heterogeneity  in  the  earthquake  source  regions  is  at  least 


partially  responsible  for  the  spectral  differences  and  that  in 


particular  the  downward  bant  lithosphere  acts  as  a waveguide  at 


some  azimuths  for  some  Rayleigh  wavelengths 


S2< 


c 


To  test  this  hypothesis,  a two-dimensional  scale  model  of 
an  assumed  island  arc  source  region  was  designed  and  con- 
structed. The  model  was  scaled  such  that  Rayleigh  energy 
was  in  the  freguency  range  20  to  200  kilohertz.  The  signals 
were  converted  to  digital  form  for  spectral  analysis.  Sources 
were  first  located  at  scaled  depths  of  0,  50  and  250  kilo- 
meters in  the  laterally  homogeneous  portion  of  the  model. 
Spectra  from  these  events  agree  guite  well  with  theory  and 
give  confidence  that  results  obtained  in  the  laterally  hetero- 
geneous portion  of  the  model,  where  theory  is  not  directly 
applicable,  are  valid.  Sources  were  then  located  at  similar 
depths  in  the  dipping  portion  of  the  model  lithosphere.  Com- 
parison of  these  spectra  with  those  from  sources  at  the  same 
depth  in  the  laterally  homogeneous  portion  indicates  a signifi- 
cant enhancement  of  energy  at  periods  less  than  50  seconds  for 
the  sources  in  the  dipping  lithosphere.  At  a period  of  30 
seconds  there  is  about  16  dB  more  amplitude  for  a source  in 
the  downward  bent  portion  of  the  lithosphere  as  compared  to 
the  same  source  at  the  same  depth  in  the  laterally  homogen- 
eous portion  of  the  model.  There  is  7 dB  more  amplitude  at 
40  seconds  and  5 dB  more  at  50  seconds.  There  is  no  differ- 
ence at  periods  of  60  seconds  and  greater.  It  is  suggested 
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that  wave  lengths  as  large  as  twice  the  width  of  the  litho 
sphere  are  utilizing  it  as  a waveguide  to  the  surface. 
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INTRODUCTION 


The  relations  for  the  far-field  surface  wave  displace- 
ments for  sources  located  in  specified  laterally  homogeneous 
Earth  models  were  presented  by  Harkrider  (1964)  and  Ben- 
Menahem  and  Harkrider  (1964),  and  have  been  expanded  upon  in 
the  work  of  Harkrider  and  Anderson  (1966)  and  Harkrider  (1970) . 
If  the  source  mechanism  and  the  elastic  parameters  of  the 
model  layers  are  specified,  the  theory  yields  smoothly-vary- 
ing Rayleigh  spectra  which  are  quite  sensitive  to  depth  of 
focus . 

There  are  data  available,  however,  which  suggest  that 
Rayleigh  wave  spectra  are  considerably  more  complicated  than 
that  predicted  by  theory.  During  the  spring  and  summer  of 
1970,  Southern  Methodist  University,  in  cooperation  with 
Teledyne-Geotech  of  Garland,  Texas,  operated  and  recorded 
digitally  a long-period  vertical  seismograph  in  a sale  mine 
near  Grand  Saline,  Texas.  The  combination  of  a broad 
response  peaking  at  45  seconds  and  the  isolation  from  atmos- 
pherically generated  noise  achieved  by  the  subterranean 
location  allowed  the  recording  of  Rayleigh  waves  of  superior 
quality.  These  seismograms  offered  a previously  unavailable 


I 


view  of  the  dispersed  wave  train  over  the  period  range  10  to 
100  seconds.  A point  transform  technique  (Filon,  1928)  was  used 
to  calculate  the  Fourier  spectra  of  the  dispersed  Rayleigh  wave 
trains  originating  from  sources  with  a wide  range  of  magnitudes, 
focal  depths,  epicentral  distances,  and  azimuths.  Many  of  these 
spectra,  such  as  the  typical  one  shown  in  Figure  1,  are  rather 
complicated  with  many  rapid  changes  in  slope--quite  different  from 
the  smoothly  varying  spectra  predicted  by  theory.  None  of  the 
Grand  Saline  spectra  is  corrected  for  instrument  response. 

A study  of  these  spectra  also  indicates  a surprising  lack  of 
attenuation  for  periods  in  the  range  20  to  50  seconds  with  in- 
creasing depth  of  focus  in  some  island  arcs. 

Figure  2 shows  the  spectra  from  two  earthquakes  with  differ- 
ent  source  depths  near  New  Britain  Island.  The  spectra  are 
normalized  at  a period  of  40  seconds.  An  increase  in  the  focal 
depth  of  from  65  kilometers  to  203  kilometers  results  in  a shift 
of  the  spectral  peak  to  longer  periods.  This  type  result  is 
also  found  for  earthquakes  in  the  New  Hebrides  and  Solomon 
Islands,  and  is  consistent  with  intuitive  expectations  of 
Rayleigh  spectral  change  with  increasing  source  depth. 

A much  different  pattern  is  seen  in  Figure  3,  which  shows 
Rayleigh  spectra  from  two  earthquakes  of  different  source  depth 


-5- 


r 

in  the  Tonga  Islands.  These  spectra  are  also  normalized  at 
a period  of  40  seconds  . In  this  case,  an  increase  in  source 
depth  from  74  kilometers  to  275  kilometers  produces  no  signi- 
ficant spectral  shift.  Figure  4 shows  the  spectra  of  Rayleigh 
waves  from  two  earthquakes  in  the  Mariana  Islands,  one  with  a 
focal  depth  of  67  kilometers  and  the  other  with  a focal  depth 
of  602  kilometers.  The  spectra  are  normalized  at  a period  of 
40  seconds.  A surprisingly  small  spectral  shift  to  longer 
periods  is  noted  for  the  deeper  focus  event.  The  Grand  Saline 
data  show  a similar  lack  of  shift  of  Rayleigh  spectra  with  in- 
creasing focal  depth  in  the  Kermedec  and  Phillipine  Island 
arcs.  No  focal  mechanism  studies  have  been  made  for  the  earth- 
quakes recorded  at  Grand  Saline,  and  thus  the  effect  that 

% 

differences  in  source  radiation  function  have  in  establishing 
this  pattern  is  not  known.  However,  it  seems  possible  that 
factors  not  included  in  layered  media  theory  are  exerting  an 
influence  in  some  source  regions.  A likely  contributing 
factor  is  lateral  heterogeneity  along  the  propagation  path 
resulting  in  (a)  horizontal  multi-pathing  and  interference, 
and  (b)  refraction,  reflection,  and  scattering  effects  in  the 
vertical  plane  containing  the  source  and  receiver.  In  this 
second  category  is  the  effect  of  the  geometry  of  the  source 
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area  with  which  this  paper  is  concerned. 

Earthquakes  are  now  known  to  originate  in  anomalous  areas 
of  the  Earth  which  are  characterized  by  lateral  heterogeneity. 

In  particular,  island  arcs  are  the  source  of  a large  majority 
of  all  earthquakes  and  the  source  of  essentially  all  deep  earth- 
quakes. Recent  evidence  from  many  fields  indicates  that  island 
arcs  are  areas  in  which  a portion  of  the  upper  mantle  is  being 
thrust  under  the  adjacent  area,  forming  a relatively  cool,  high- 
strength  zone  in  the  midst  of  hotter,  weaker  material.  Earth- 
quakes originate  along  this  downward-bent  wedge  of  lithosphere, 
suggesting  that  a more  realistic  model  for  surface  wave  spectral 
studies  is  that  of  a source  at  some  depth  below  an  elbow  in  a 
curved  waveguide. 

It  is  interesting  to  note  that  the  Marianas,  Tonga,  Kermedec, 
and  Phillipine  Islands,  where  earthquakes  recorded  at  Grand  Saline 
show  no  large  shift  of  Rayleigh  spectral  peak  to  longer  periods  with 
increasing  depth  of  focus,  are  subduction  zones  in  which  the  litho- 
sphere plunges  away  from  the  recording  station  at  Grand  Saline.  On 
the  other  hand,  the  New  Britain,  New  Hebrides,  and  Solomon 
Islands,  where  earthquakes  do  show  a spectral  shift,  are  zones 
in  which  the  lithosphere  plunges  toward  Grand  Saline.  This 
suggests  that  the  geometry  of  the  dipping  lithosphere  may  be 
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enhancing  the  Rayleigh  periods  of  20  to  50  seconds  at  those 
azimuths  opposite  to  the  direction  of  dip. 

An  Earth  model  in  which  a downgoing  slab  of  lithosphere 
is  considered  greatly  complicates  the  theoretical  problem. 
Numerical  solutions  to  the  laterally  heterogeneous  problem 
are  possible  in  principle,  but  would  be  complicated  and  time 
consuming  in  practice.  As  an  alternative  approach,  a two- 
dimensional  analog  model  of  an  assumed  island  arc  structure 
with  a downgoing  wedge  was  designed  and  constructed,  and  the 
effect  of  this  source  structure  on  the  Rayleigh  spectrum  was 
evaluated  using  the  scale  model. 

Several  studies  have  been  made  previously  of  the  propa- 
gation of  Rayleigh  waves  in  two-dimensional  models.  These 
have  been  primarily  nor  scale  models  of  Earth  situations,  but 
rather  observations  of  the  behavior  of  Rayleigh  waves  as  they 
propagate  around  corners,  across  non-horizontal  boundaries, 
and  through  channels  of  changing  width.  For  example,  de 
Bremaecker  (1950)  studied  the  energy  conversions  occurring 
when  a Rayleigh  wave  is  incident  upon  a corner  whose  angle 
is  between  0 and  180  degrees.  Knopoff  and  Gangi  (1961),  in 
a similar  experiment,  observed  that  the  Rayleigh  wave  under- 
goes a change  in  wave  form  as  it  negotiates  the  corner.  Kuo 
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and  Thompson  (1963)  studied  experimentally  the  effect  of  a 
gently-sloping  interface  on  the  propagation  of  Rayleigh  waves. 
Gangi  (1967)  determined  experimentally  the  P wave  to  Rayleigh 
wave  conversion  coefficients  at  a stress-free  wedge.  Ottaviani 
(1971)  showed  that  P waves  incident  upon  a corner  generate 
secondary  Rayleigh  waves  which  then  propagate  in  either  direc- 
tion from  the  corner.  All  of  these  investigators  used  equip- 
ment and  procedures  based  on  a work  by  Oliver,  Press,  and 
Ewing  (1954)  who  demonstrated  mathematically  and  verified 
empirically  that  wave  propagation  in  ana  on  the  edge  of  thin 
sheets  can  be  related  to  plane  wave  propagation  in  three  de- 
mens ions.  The  characteristic  relation  between  period  and 
phase  velocity  for  Rayleigh  waves  on  the  edge  of  a plate  is 
the  same  as  for  the  free  surface  of  an  infinite  solid.  How- 


ever, the  plate  dilatational  velocity  is  Vp 


rather  than  the  infinite  solid  dilatational  velocity, 

. I i\  I JLA.  I / i 

ek.  — 


■AJhlfcLT 

. P J 


, where  A andyi.  are  the  Lame  elastic  con- 
stants and  I?  is  the  density.  By  simply  replacing  J*  = Vp, 
many  of  the  problems  of  propagation  of  plane  waves  in  a 


stratified  Earth  can  be  related  to  two-dimensional  models 
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CONSTRUCTION  OF  THE  MODEL 


In  trying  to  build  a scale  model  of  a portion  of  the 
Earth  it  is  quickly  apparent  that  the  classical  technique  of 
two-dimensional  model  construction,  i.e.,  butt-joining  pieces 
of  different  materials  to  fashion  layers  of  different  elastic 
properties,  is  not  satisfactory.  There  is  not  a wjde  enough 
assortment  of  available  modeling  materials  to  build  even  the 
grossest  approximation  of  a crustal  and  upper  mantle  structure. 
Table  1 shows  most  of  the  commonly  available  modeling  materials 
and  their  elastic  velocities.  While  it  is  easy  enough  to  find 
three  or  four  materials  with  significantly  different  seismic 
velocities,  it  is  quite  another  matter  to  obtain  three  or  four 
materials  whose  velocities  differ  by  a desired  specific  ratio, 
which  is  the  necessary  criterion  for  building  a scale  model. 

A solution  to  this  problem  was  suggested  by  Angona  (1960) 
and  Riznichenko,  et  al . (1961).  If  sheets  of  two  materials 

are  joined  broadside  to  form  a composite  plate,  the  seismic 
velocities  of  wave  lengths  long  in  comparison  to  the  total 
thickness  will  be  intermediate  between  the  velocities  of  the 
component  plates  and  will  be  proportional  to  the  relative 
thicknesses  of  the  components.  The  exact  empirical  relation 


TABLE 


-i  1 

IV  l,f><  + <-ipi 


where  f\y  is  the  composite  velocity 

is  the  thickness  proportion  of  com- 
ponent i 

is  the  velocity  of  component  i 
p i is  the  density  of  component  i 
By  selecting  a high  velocity  material  such  as  aluminum, 
and  low  velocity  material  such  as  vinyl,  a wide  range  of  velo- 
cities can  be  obtained  by  varying  the  relative  thicknesses. 

The  choice  of  the  particular  structure  of  the  crust  and 
upper  mantle  to  be  modeled  was  guided  by  a desire  to  have  a 
suitably  Earth-like  situation  incorporating  first-order  features 
while  avoiding  unnecessarily  detailed  layering.  Thus,  only  four 
Earth  units  were  modeled — a crust,  a high-velocity  lithosphere 
with  a dipping  wedge,  a low-velocity  intermediate  layer,  and  a 
high-velocity  half-space.  The  Earth  shear  wave  velocity  dis- 
tribution, shown  in  Figure  5,  which  is  a simplification  of 
the  oceanic  shear  wave  model  of  Toksoz  and  Anderson  (1966), 
was  adopted  for  modeling. 

Aluminum  has  the  highest  elastic  velocities  of  any  easily 
available  modeling  material.  Tne  basic  scaling  problem  in 
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constructing  the  model  was  to  utilize  the  aluminum  as  the 
high-velocity  half-space,  together  with  varying  thicknesses 
of  vinyl  to  form  the  lower  velocity  layers.  Since  the  shear 
wave  velocity  of  aluminum  is  3.24  km/sec  and  the  desired  half- 
space velocity  is  5.4  km/sec,  a velocity  scale  factor  of  5/3 
was  suggested.  That  is 

or  Earth  = 5/3  nr 

Model 

The  frequency  range  of  interest,  (0.01  to  0.1  Hz),  can 
be  scaled  to  the  effective  frequency  range  of  the  piezoelectric 
transducers  (20  to  200  kHz).  This  implies  a time  (t)  scale 
factor  such  that 


L Earth  = 2 x 10^  L Model 

By  fixing  the  time  and  velocity  scale  factors,  the  length  re- 
lation is  constrained  to  be 


/ Earth  = Earth  • ^ Earth  = (10/3  x 106)  s-  Model 

The  relative  thicknesses  of  the  vinyl  and  the  aluminum  in 
the  composite  sheets  forming  the  upper  zones  were  determined 
by  equation  1.  A diagram  showing  the  dimensions  of  the  model 
and  the  widths  and  velocities  of  each  layer  are  shown  in 
Figure  6.  Shown  in  Figure  7 is  a schematic  cross-sectional 
view  of  the  composite  plate  model.  The  four  different  velocity 
layers  were  made  from  different  combinations  of  two  thicknesses 
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of  vinyl,  0.006  inches  and  0.016  inches.  Each  layer  of  vinyl 
was  attached  with  Duro  E*Pox*E  resin,  number  14. 

The  maximum  thickness  of  the  composite  plate  was  0.104  in- 
ches, or  equivalently,  0.264  centimeters.  The  shortest  wave 
length  of  interest  was  2.39  centimeters.  It  is  important  that 
the  thickness  of  a two-dimensional  model  be  small  compared  to 
the  shortest  wave  lengths  of  interest.  Otherwise,  the  com- 
press ional  wave  will  be  dispersive,  which  of  course  is  not 
analogous  to  compress ional  wave  propagation  in  the  Earth. 

The  Source 

A piezoelectric  crystal  was  used  as  the  source  of  the 
elastic  waves.  This  type  of  crystal  changes  shape  when  an 
electric  field  is  applied,  and  conversely,  develops  an  electric 
field  when  pressure  is  applied.  Rochelle  salt,  ammonium  dehy  — 
drogen  phosphate,  lithum  sulphate,  barium  titanate  crystals, 
and  lead  z irconate-lead  titanate  (PZT)  ceramics  are  in  common 
use  as  electromechanical  transducers. 

Various  types  of  PZT  ceramics  are  available,  including 
radial  expanders,  length  expanders,  shear  actors,  and  bimorphs 
(benders) . Each  type  was  investigated  for  the  purpose  of  de- 
termining which  was  most  suitable  for  the  purposes  of  the  study. 
The  chief  criteria  involved  in  the  selection  and  implementation 
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of  the  source  were  that  it  must  be  a type  which  could  effi- 
ciently gopa-ate  surface  waves  when  located  at  depth  Waa 
attachable  to  the  plate  in  some  fashion  so  as  to  excite 
symmetrical  modes,  i.e.,  to  induce  energy  along  the  medial 
plane  of  the  plate  and  not  transverse  to  it.  Length  expanders 
were  found  to  be  satisfactory  in  every  respect.  The  type  used 
was  1/8-inch  long  and  acted  as  a vector  force  directed  along 
the  line  of  orientation  of  the  crystal.  The  crystal  was 
mounted  in  a 3/16-inch  diameter  hole  drilled  through  the  plate. 
The  crystal  was  held  against  the  interior  of  the  hole  by  a 
rubber  band  which  was  stretched  prior  to  insertion  of  the 
crystal  and  then  allowed  to  contract,  the  increased  thickness 
holding  the  crystal  firmly. 

The  1/8-inch  Clevite  PZT-5  piezoelectric  crystal  has  a 
resonant  frequency  near  350  kilohertz.  If  this  frequency  falls 
within  the  passband  of  interest,  some  measure  must  be  taken  to 
dampen  the  resonant  ringing.  This  can  be  done  by  backing  the 
crystal  with  ordinary  modeling  clay.  In  the  present  study, 
the  resonant  frequency  was  well  outside  the  20  to  200  kilohertz 
passband  and  was  not  evident  on  the  seismograms;  thus  the  bulky 
backing  was  not  necessary  and  only  the  1/8-inch  crystal  itself 
needed  to  be  applied  to  the  interior  of  the  small  hole. 


In  implanting  the  crystal  in  the  plate,  some  care  was 


, 

. 
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taken  to  place  the  crystal  face  squarely  against  the  hole 
edge  and  in  good  contact  with  the  aluminum  portion  of  the 
composite  plate.  Otherwise,  asymmetric  or  low-level  signals 
could  be  generated.  Good  contact  was  aided  by  filing  slightly 
the  interior  of  the  hole  to  remove  the  curvature  at  the  point 
of  contact. 

Since  making  a hole  in  the  model  was  necessary  to  implant 
a source  at  depth,  a preliminary  experiment  was  conducted  to 
determine  the  effect  of  a previously  drilled  hole  in  the 
vicinity  of  a source  location.  A 3/16-inch  diameter  hole  was 
drilled  in  a 48  inches  by  36  inches  by  0.040  inch  sheet  of 
aluminum,  the  center  of  the  hole  being  1.0  inch  from  one  edge. 
A horizontal  source  was  mounted  in  the  hole,  and  a seismogram 
was  recorded.  A second  hole  was  then  drilled  2.25  inches  to 
the  side  of  the  first  hole,  but  not  in  the  travel  path.  The 
source  was  then  located  in  the  second  hole  and  another  seismo- 
gram was  recorded.  The  resulting  seismograms  were  essentially 
identical,  indicating  that  any  reflected  or  scattered  energy 
resulting  from  the  superfluous  hole  is  not  measurable  at  the 
operating  system  amplification. 
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It  was  the  procedure  to  conduct  the  main  body  of  experi- 
ments in  a sequence  such  that  no  hole  was  ever  located  along 
a travel  path,  but  was  always  located  to  the  rear  of  the  source, 
as  in  the  described  preliminary  experiment. 

The  Receiver 

A Clevite  1/8-inch  PZT-4  ceramic  crystal  was  used  as  the 
sensing  element  and  was  backed  and  amplified  as  shown  in  Fig- 
ure 8.  The  crystal  and  amplifier  were  mounted  in  a small 
carriage  which  rode  upon  the  upper  edge  of  the  model.  Gravity 
was  sufficient  to  hold  the  crystal  in  good  contact  with  the 
plate  edge,  although  care  was  taken  to  place  the  crystal  face 
squarely  against  the  plate  edge.  Failure  to  do  this  will  re- 
sult in  the  loss  of  some  high  frequency  energy. 

The  Modeling  Electronics 

A block  diagram  of  the  modeling  electronics  is  shown  in 
Figure  9.  The  piezoelectric  crystal  is  activated  by  a 
Thyratron  pulser  which  generates  a step— function  voltage  at  a 
prescribed  (but  variable,  if  desired)  rate.  The  positive 
voltage  step  causes  the  ceramic  transducer  to  expand  and  thus 
to  generate  elastic  waves  at  the  point  of  contact  with  the 
plate.  The  Thyratron  is  set  to  generate  a new  signal  at  a 
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rate  of  several  tens  of  times  per  second.  A delayed  version 
of  the  Thyratron  excitation  pulse  is  used  to  trigger  a sampling 
oscilloscope  which  displays  the  signal  after  recording.  Prior 
to  oscilloscope  display,  the  recorded  signal  is  amplified  and 
filtered  so  that  the  effective  passband  is  20  to  200  kilohertz. 
Thus,  signal  periods  are  in  the  range  of  5 to  50  microseconds, 
and  travel  times  are  generally  a few  hundred  microseconds.  A 
paper  copy  of  the  seismogram  on  an  expanded  time  scale  was  rou- 
tinely obtained  by  means  of  a wave  translation  feature.  Upon 
activation  of  a switch,  a sampling  pulse  proceeds  to  sweep  across 
the  stored  signal  at  a slow,  regulated  speed  which  is  coordinated 
with  the  paper  chart  speed.  The  resulting  paper  seismogram  has 
a time  base  of  42.86  /Js  per  inch. 

The  influence  of  the  frequency  responses  of  the  source, 
the  receiver,  and  of  the  filtering  and  amplification  electronics 
on  the  spectrum  of  the  model  Rayleigh  signal  can  be  seen  in 
Figure  10.  This  spectrum  was  obtained  by  placing  the  source 
and  the  receiver  immediately  adjacent  to  each  other,  thus  eli- 
minating the  response  of  the  medium  from  the  system.  Also  shown 
in  Figure  10  is  the  spectrum  from  a typical  seismogram  recorded 
at  Grand  Saline.  The  frequency  response  of  the  seismograph 
has  not  been  removed.  Note  that  the  model  spectrum  appears  to 
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be  a close  approximation  to  a smoothed  version  of  the  actual 
spectrum.  This  indicates  that  the  model  transducers  and 
electronics  produced  a spectrum,  which  when  adjusted  by  the 
model  time  scale  factor,  encompasses  the  frequency  range  of 
interest . 

Digitization 

The  wave  transformed  output  was  routed  to  a digitizer 
which  war ii  terfaced  with  an  XDS-925  computer.  The  analog  data 
were  amplified  by  a factor  of  10  to  boost  the  signal  into  the 
+ 10  volt  range  required  by  the  digitizer,  and  were  passed 
through  an  anti-aliasing  filter  to  remove  noise  at  a scaled 
period  below  10  seconds.  A sampling  trigger  was  supplied  to 
the  digitizer  from  the  Thyratron  pulser  after  delaying  and  con- 
ditioning the  pulse.  A 1200  /Os  delay  was  introduced  to  allow 
for  settling  time  of  the  sample-and-hold  amplifier.  The  delayed 
pulse  was  conditioned  to  provide  a sampling  trigger  of  +10  volts 
amplitude  and  2.5  JUs  duration. 

By  suitably  adjusting  the  Thyratron  repitition  rate  and 
allowing  for  the  time  scale  factor  and  the  wave  translation 
time  factor,  the  resultant  signal  was  sampled  once  per  second. 

Spectral  Determinations 


The  digitized  seismograms  were  written  on  magnetic  tape 


in  a format  compatible  with  data  processing  programs  which 
were  already  being  used  on  the  real  seismograms  recorded  at 
Grand  Saline.  One  of  these  programs,  a point  transform  tech- 
nique (Filon,  1928)  was  used  to  calculate  the  Rayleigh  spectra 
obtained  in  this  study.  This  method  evaluates  the  Fourier 
integral  of  the  desired  time  series  by  dividing  the  time  series 
into  2N  equal  parts  with  an  interval  h equal  to  the  sampling 
interval.  It  is  assumed  that  over  the  2N  range  the  time  series 
can  be  approximated  with  sufficient  accuracy  by  a parabola. 
Since  the  parabola  is  constrained  at  three  points  (the  three 
sample  values),  the  arbitrary  constants  or  the  parabola  are 
uniquely  determined  and  can  be  expressed  in  terms  of  the  adja- 
cent sample  values  and  the  sample  interval.  ine  integrals  over 
each  2N  range  can  then  be  evaluated  analytically  by  parts  and 
can  be  combined  to  give  the  total  integral.  Contrary  to  the 
Fast  Fourier  Transform  (e.g.,  Welch,  1967),  the  Filon  tech- 
nique allows  the  transform  to  be  calculated  at  any  desired 
frequency.  No  smoothing  function  was  applied  to  the  spectra 
since  the  Rayleigh  wave  is  assumed  to  be  a unique  entity,  hav- 
ing a beginning  and  an  end,  and  not  to  be  a sample  selected 
from  an  infinite  set  as  is  usually  the  assumption  in  spectral 
analyses  of  seismic  time  series. 


After  computation,  the  spectra  were  automatically  plotted 
by  a 7126  XDS  plotter.  The  specific  quantity  plotted  was  the 
square  root  of  the  sum  of  the  squares  of  the  real  and  the 
imaginary  parts  of  the  transform  expressed  in  decibels  rela- 
tive to  a convenient  value  fixed  by  the  dynamic  range  of  the 
digital  recording  system. 

Preliminary  Experiments 

The  aluminum  alloy  used  to  build  the  model  was  2024-T3. 

The  Rayleigh  velocity  of  the  alloy  was  determined  by  measuring 
the  travel  times  over  several  distances.  The  reciprocal  of  the 
slope  of  the  best  fit  straight  line  was  taken  as  the  Rayleigh 
velocity.  However,  the  scatter  of  the  data  points  was  so  small 
that  subsequent  velocity  measurements  were  made  at  only  one 
distance . 

To  verify  formula  1,  composite  plates  of  aluminum  and 
vinyl  were  constructed  for  each  of  the  three  thickness  combina- 
tions shown  in  Figure  7.  The  elastic  velocities  shown 
in  Figure  6 were  then  determined  by  measuring  the  travel 
times  over  a 400  millimeter  distance. 

A surface  source  on  an  aluminum-vinyl  half-space  dupli- 
cates the  classical  theoretical  problem  of  Lamb  (1904) . Figure 
11  shows  the  model  seismogram  and  Lamb's  theoretical  result. 
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The  close  correspondence  of  the  two  indicates  that  the  composite 
plate  does  indeed  set  rr  h V»r1 £ = ppRtia . Tnat:  ip,  4-r  no 

modal  propagation  within  the  separate  vinyl  layers,  no  body 
wave  dispersion  which  would  result  if  the  total  plate  thickness 
were  large  compared  to  the  smallest  wave  lengths,  nor  is  there 
evidence  of  asymmetrical  mode  propagation.  To  demonstrate  the 
undesirable  asymmetric  mode,  the  thickness  expander  source  was 
applied  transversely  to  the  aluminum  plate.  Figure  12  shows 
the  result.  The  very  long  period  inversely  dispersed  wave  train 
is  unmistakable. 

It  is  interesting  to  compare  the  model  results  for  a source 
at  depth  within  a homogeneous  half-space,  represented  by  an 
aluminum  plate,  with  Lapwood's  (1949)  theoretical  results.  Lap- 
wood  considered  a cylindrical  explosion  as  a source  at  a depth 
h.  Figure  13  shows  a model  seismogram  for  the  case  where  both 
shear  and  compress ional  waves  are  generated  by  the  source  and 
have  travelled  a distance  x which  is  large  compared  to  depth  h. 
The  fact  that  the  source  is  located  below  the  surface  intro- 
duces several  phase  arrivals  which  are  not  present  in  Lamb's 
problem.  The  phase  denoted  as  sP  has  travelled  from  the  source 
to  the  surface  as  a shear  wave  and  has  then  been  refracted  a- 
long  the  surface  as  a compressional  wave.  The  phase  pS  travelled 
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to  the  surface  as  a compressional  wave  and  along  the  surface 
as  a shear  wave.  The  arrival  times  of  each  phase  as  given 
by  Lapwood  in  terms  of  x(  h,  and  the  shear  and  compressional 
velocities  of  the  half-space  are  indicated  in  the  figure. 

Each  of  the  theoretically  predicted  phases  can  be  observed 
on  the  model  seismogram  and  are  seen  to  be  arriving  at  the 
times  determined  by  Lapwood.  Since  the  theoretical  source 
space  and  time  function  are  different  from  the  model  source, 
no  exact  similarity  can  be  expected  between  the  corresponding 
phase  wave  forms. 


RAYLEIGH  SPECTRA  ORIGINATING  FROM  FOCI  IN 
LATERALLY  HOMOGENEOUS  LAYERED  MEDIA 

As  can  be  seen  from  the  sketch  of  the  model  in  Figure  6, 
it  was  possible  to  place  sources  in  the  laterally  homogeneous 
portion  of  the  model,  as  well  as  in  the  downward  bent  portion 
of  the  model  lithosphere.  This  was  done  in  order  to  establish 
a set  of  standard  spectra  with  which  to  compare  the  spectra 
originating  from  the  lithospheric  wedge,  and  also  to  compare  in 
a qualitative  way  with  theoretical  spectra. 

Rayleia";  waves  were  recorded  at  a distance  of  800  milli- 
meters, aled  distance  of  2667  kilometers,  from  sources  at 

the  surf,  e,  and  at  scaled  depths  of  50  and  250  kilometers.  The 
distance  of  800  millimeters  was  convenient  for  two  reasons.  It 


was  short  enough  that  no  compress ional  wave  reflections  from  the 
plate  edges  could  interfere  with  the  Rayleigh  wave  arrival . It 
was,  however,  several  times  larger  than  the  longest  wave  lengths 
of  interest  and  therefore  allowed  comparison  of  the  results  with 
far-field  theory. 

Figure  14  shows  the  seismogram  of  a Rayleigh  wave  from  a 
vertical  surface  force.  The  spectrum  of  the  Rayleigh  wave  is 
shown  in  Figure  15.  The  time  window  from  which  the  signal 


spectrum  was  taken  is  shown  in  Figure  14.  Above  the  noise  the 
spectrum  is  perfectly  smooth,  but  at  a level  where  the  noise 
dominates,  the  spectrum  becomes  unstable.  Note  that  the  signal- 
to-noise  ratio  is  about  30  dB  near  the  middle  of  the  pass-band. 

It  is  interesting  to  compare  the  Rayleigh  spectrum  with  a 
theoretically  predicted  spectrum  for  a laterally  homogeneous 
oceanic  structure.  Such  a spectrum  is  presented  by  Harkrider 
and  Anderson  (1966)  for  a vertical  surface  source.  For  compara- 
tive purposes  the  theoretical  spectrum  must  be  adjusted  for  the 
model  source  spectrum  and  for  the  model  attenuation  as  a func- 
tion of  frequency  and  distance.  By  comparing  the  spectrum  in 
Figure  15  with  the  source  spectrum  in  Figure  10  the  attenuation 
with  distance  of  each  frequency  can  be  determined.  The  spectra 
are  the  same  at  periods  greater  than  55  seconds,  but  at  smaller 
periods  there  is  significant  attenuation  with  travel  path.  Ta- 
ble 2 shows  amplitude  attenuation  in  dB  as  a function  of  scaled 
travel  path  distance.  By  applying  these  attenuation  factors  to 
the  Harkrider  and  Anderson  theoretical  spectrum  and  adjusting 
the  resultant  for  the  model  source  spectrum,  the  theoretical 

spectra  can  be  directly  compared  to  the  model  spectrum  for  the 
same  depth. 

The  dotted  curve  in  Figure  15  shows  the  theoretical  spectrum 
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Period 

(sec) 

Attenuation 

(dB/km) 

20 

. 01900 

30 

.00600 

40 

. 00075 

50 

. 00033 

55 

. 00000 

60 

. 00000 

70 

.00000 

80 

. 00000 

90 

. 00000 

100 

.00000 
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for  the  surface  vertical  force.  Figures  16,  17  and  18  show 
similar  comparisons  for  a horizontal  force  at  a depth  of  50 
kilometers,  a vertical  force  at  50  kilometers,  and  a horizon- 
tal force  at  250  kilometers,  respectively.  The  spectra  are 
normalized  at  a period  of  40  seconds.  The  spectra  shown  in 
Figure  16  are  particularly  interesting.  The  theoretical  spectrum 
for  the  horizontal  force  at  a depth  of  50  kilometers  shows  a zero 
at  a period  of  60  seconds.  The  model  spectrum  also  shows  the  zero, 
but  it  has  been  filled  in  to  the  noise  level;  the  noise  also 
causes  a slight  shift  in  the  period  at  which  the  'hole'  appears 
to  occur.  The  secondary  peaks  of  the  two  spectra,  occurring  at 
a period  of  about  85  seconds,  match  very  well.  It  should  be  noted 
that  the  theoretical  spectra  are  computed  on  the  basis  of  a 
very  detailed  version  of  the  Toksoz  and  Anderson  (1966)  oceanic 
structure  while  the  seismic  model  is  merely  a four-layer  simpli- 
fication. Thus,  the  comparison  of  theoretical  and  model  spectra 
cannot  be  considered  as  a rigorous  demonstration  of  the  fidelity  of 
the  model  construction.  Nevertheless,  the  very  close  agreement 
between  the  model  spectra  and  the  theoretical  spectra  for  all 
source  depths  indicates  that  surface  wave  generation  and  propaga- 
tion  in  the  laterally  homogeneous  portion  of  the  model  are  correctly 
scaled  to  represent  phenomena  in  an  oceanic  structure. 
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SPECTRA  OF  RAYLEIGH  WAVES  FROM  SOURCES 
IN  THE  LATERALLY  HETEROGENEOUS  ZONE 

Sources  were  located  at  scaled  depths  of  0,  50,  150,  and 
250  kilometers  in  the  downward  bent  portion  of  the  model  litho- 
sphere. Figure  19  shows  the  source  locations  in  relations  to 
the  geometry  of  the  anomalous  zone.  Note  that  the  line  of  foci 
is  somewhat  to  the  top  of  center  of  the  lithosphere.  It  has 
been  proposed  that  all  earthquakes  occur  in  the  upper  20  kilo- 
meters (Isacks,  Oliver,  and  Sykes,  1968)  of  the  downgoing  slab, 
and  alternatively  that  the  earthquakes  occur  nearer  the  center 
a ‘■ong  the  a*is  of  cooler  temperatures  (e.g.,  McKenzie,  1969;  Griggs, 
1972) . There  does  not  appear  to  be  sufficient  evidence  to  dis- 
tinguish between  the  two  suggestions  at  this  time.  In  this 
study  the  latter  hypothesis  was  followed,  although  it  is  un- 
likely that  it  would  make  much  difference  since  the  spectra  for 
sources  at  0 and  50  kilometers  depth  in  the  laterally  homogene- 
ous zone  are  almost  the  same. 

A Rayleigh  spectrum  from  a vertical  surface  force  at  a 
scaled  distance  of  2667  kilometers  is  shown  in  Figure  20.  There 
appears  to  be  no  difference  between  this  spectrum  and  the 
spectrum  of  the  surface  source  in  the  laterally  homogeneous 
portion  of  the  model,  which  is  shown  in  the  same  figure  for 
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comparison.  Similarly,  the  Rayleigh  spectrum  from  a vertical 
force  above  the  dipping  lithosphere  at  a scaled  depth  of  50 
kilometers,  shown  in  Figure  21  is  very  similar  to  the  analogous 
spectrum  from  a;  source  at  a depth  of  50  kilometers  in  the  later- 
ally homogeneous  area,  which  is  also  shown  in  the  figure. 

Figure  22  shows  the  spectrum  of  a Rayleigh  wave  from  a 
horizontal  force  at  a scaled  depth  of  250  kilometers  in  the  ano- 
malous zone.  Also  shown,  by  alternating  dots  and  dashes,  is  the 
theoretical  curve  for  the  same  source  in  laterally  homogeneous 
media  and,  by  a solid  line,  the  model  results  for  the  same  source 
in  laterally  homogeneous  media.  The  solid  and  the  dot-dash  spectra 
are  almost  coincident  in  the  passband  30  to  100  seconds,  indicat- 
ing that  the  model  and  theoretical  results  for  laterally  homo- 
geneous source  areas  and  travel  paths  are  in  agreement.  Below 
a period  of  30  seconds,  the  model  signal  spectrum  becomes  lost 
in  the  noise.  Comparison  of  the  solid-lined  spectrum  of  the 
source  in  the  downward  bent  lithosphere  to  the  other  two  spectra 
indicates  a significant  enhancement  of  energy  in  the  period 
range  30  to  60  seconds.  With  the  amplitudes  equcilized  along 
the  long-period  side  of  the  spectra,  there  is  about  16  dB  more 
amplitude  for  the  dipping  lithosphere  source  at  a period  of  30 
seconds.  There  is  7 dB  more  amplitude  at  40  seconds  and  5 dB 


more  at  50  seconds,  with  the  curves  converging  at  60  seconds. 

The  source  depth  of  250  kilometers  is  well  below  the  bend  in 
the  lithosphere,  and  it  appears  that  wave  lengths  as  large  as 
twice  the  width  of  the  lithosphere  are  utilizing  it  as  a wave- 
guide to  the  surface. 

Focal  mechanism  studies  of  earthquakes  in  island  arcs  (Isacks, 
Oliver,  and  Sykes,  1968)  indicate  that  the  compressive  axes  of 
earthquakes  occurring  below  the  bend  in  the  lithosphere  are  gen- 
erally aligned  along  the  dip  of  the  lithosphere.  Accordingly,  a 
piezoelectric  crystal  was  situated  at  a scaled  depth  of  250  kilo- 
meters so  that  the  axis  of  compression  was  along  the  dip  of  the 
model  lithosphere.  The  resulting  Rayleigh  spectrum  is  shown  in 
Figure  23,  where  it  is  compared  with  the  250  kilometer  source  in 
the  layered  half-space.  This  orientation  of  the  source  also  shows 
the  spectral  shift  to  higher  frequencies  caused  by  the  location  of 
the  source  in  the  anomalous  zone  relative  to  the  location  in  the 
layered  half-space. 

To  determine  the  behavior  of  the  Rayleigh  spectrum  nearer  to 
the  bend  in  the  lithosphere,  a source  was  situated  at  a scaled 
depth  of  150  kilometers  in  the  anomalous  zone.  The  Rayleigh 
spectrum  is  shown  in  Figure  24,  along  with  previously  shown 
spectra  for  sources  at  50  kilometers  and  250  kilometers  deep  in 
the  anomalous  zone.  The  bend  in  the  model  lithosphere  occurs 
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ci  t a depth  of  100  kilometers,  so  that  Figure  24  represents 
spectra  from  sources  50  kilometers  above  the  bend,  50  kilometers 
below  the  bend,  and  150  kilometers  below  the  bend.  This  figure 
illustrates  the  change  in  the  Rayleigh  spectrum  with  increasing 
source  depth  in  the  dipping  lithosphere  and  can  be  contrasted 
to  Figure  25  which  shows  the  spectral  shift  occurring  with  in- 
creasing source  depth  in  the  laterally  homogeneous  portion  of 
the  model.  The  spectra  shown  in  Figures  20  through  25  are  normal- 
ized at  a period  of  55  seconds.  In  the  layered  half— space  there 
is  a significant  and  theoretically  predictable  spectral  shift  to 
longer  periods  with  increasing  source  depth.  For  sources  in  the 
dipping  lithosphere,  there  is  no  corresponding  shift  of  the  spec- 
tral peak;  indeed,  the  shape  of  the  spectra  are  very  similar  at 
periods  less  than  60  seconds,  with  the  most  obvious  effect  of 
increasing  source  depth  being  the  relative  increase  in  energy  in 
the  period  range  60  to  100  seconds.  Thus,  in  the  dipping  litho- 
sphere, there  is  a relative  broadening  of  the  spectral  peak  to 
the  long  period  side  with  increasing  source  depth,  but  little 
accompanying  decrease  on  the  short  period  side,  at  least  to  a 
depth  of  250  kilometers.  The  spectral  change  with  increasing 
source  depth  illustrated  in  Figure  24  is  interpreted  to  mean 
that  wave  lengths  as  large  as  twice  the  width  of  the  lithosphere 
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are  utilizing  the  dipping  lithosphere  as  a waveguide  to  the 
surface,  while  larger  wave  lengths  are  not  enhanced  and  tend 
to  behave  as  if  the  sources  were  at  similar  depths  in  a layered 
half-space. 
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DISCUSSION 


Simple  one-directional  forces  were  used  as  the  model  sources 
to  produce  the  Rayleigh  spectra.  This  was  done  to  keep  the  source 
as  simple  as  possible  while  studying  the  effect  of  the  source  re- 
gion. However,  focal  mechanism  studies  (Isacks,  Oliver,  and 
Sykes,  1968)  show  that  a double  couple,  or  shear  dislocation,  is 
a more  realistic  source  function.  For  that  reason  a quantitative 
comparison  of  the  model  Rayleigh  spectra  with  spectra  of  real 
earthquakes  recorded  at  Grand  Saline  is  not  warranted.  Neverthe- 
less, some  interesting  qualitative  observations  can  be  made. 

As  noted  earlier,  the  Grand  Saline  spectra  are  characterized 
by  rapid  changes  of  slope  and  "holes"  within  the  signal  passband 
of  20  to  100  seconds.  The  holes  are  often  of  sufficient  relief 
to  effectively  mask  possible  depth-dependent  spectral  zeroes  re- 
sulting from  source  displacement  nodes.  The  model  spectra  are 
devoid  of  holes  other  than  the  theoretically-predicted  energy 
zero  at  60  seconds  for  a horizontal  source  at  a depth  of  50  kilo- 
meters. Nc  smoothing  function  was  used  on  either  the  real  or 
the  model  spectra.  The  model  results  indicate  that  the  holes 
in  the  real  spectra  are  not  due  to  interference  phenomena 
associated  with  source  location  in  the  anomalous  zone.  It  seems 
likely  that  the  holes  result  from  interference  effects  due  to 


lateral  refraction  along  the  propagation  path.  The  presence 
of  the  holes  could  prove  a serious  obstacle  to  the  use  of  fre- 
quency-dependent spectral  zeroes  to  determine  focal  depth. 

The  fact  that  Rayleigh  spectra  from  foci  above  the  dipping 
portion  of  the  lithosphere  are  not  significantly  different  from 
spectra  from  foci  at  the  same  depths  in  the  layered  half-space 
implies  that  current  theories  of  spectral  determination  based 
on  a laterally  homogeneous  layered  media  are  adequate  for  shallow 
events  in  island  arc  regions,  the  focal  depth  limit  of  applica- 
bility probably  being  equal  to  the  thickness  of  the  lithosphere. 
Since  the  lithosphere  thickness  may  vary  from  arc  to  arc,  the 
depth  limit  of  applicability  would  vary  accordingly. 

Because  of  size  limitations  of  the  model,  no  seismograms 
could  be  recorded  on  the  backward  side  of  the  island  arc.  The 
Grand  Saline  data  suggest  chat  the  waveguide  effect  is  not  oper- 
ative in  that  direction,  although  the  degree  of  conformity  of 
such  spectra  with  layered  media  theory  is  problematical.  At 
any  rate,  spectra  for  earthquakes  occurring  below  the  bend  in 
the  lithosphere  are  no  doubt  dependent  on  azimuth.  Surface 
wave  magnitudes  (Ms)  are  known  to  have  an  azimuthal  dependence, 
but  this  has  been  attributed  completely  to  the  radiation  pat- 
tern of  the  source  function.  While  source  radiation  pattern 


may  be  the  controlling  factor  for  shallow  earthquakes,  Lhe 
effect  of  the  geometry  of  the  source  region  must  be  considered 
for  intermediate  and  deep  focus  earthquakes.  For  a given  source 
function,  the  higher  values  of  Ms  would  be  observed  at  azimuths 
opposite  to  the  direction  of  dip  of  the  lithospheric  plate. 

The  technique  of  utilizing  a two-dimensional  model  made 
by  face  joining  two  or  more  materials  of  different  elastic 
properties  has  proven  to  be  a practical  and  flexible  method  of 
formulating  structures  of  geophysical  interest.  It  should  find 
other  applications  in  situations  where  it  is  desired  to  investi- 
gate the  effect  on  seismic  phenomena  of  zones  of  great  lateral 
structural  variation,  such  as  mid-ocean  ridges  and  continental 
boundaries . 
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